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This  technical  report  consi3ts  of  three  parts.  The  central  problem 
is  the  estimation  of  three-dimensional  motion  parameters  of  a  rigid  planar 
patch  from  image  sequences  (each  frame  is  a  central  projection). 

In  Part  I,  we  show  that  given  two  Image  frames  one  can  determine 
uniquely  (by  solving  linear  equations)  eight  "pure  parameters"  which  are 
nonlinear  functions  of  the  actual  motion  parameters.  In  Part  II,  a  method 
is  presented  for  determining  the  motion  parameters  from  the  eight  pure  para¬ 
meters.  The  method  requires  the  singular  value  decomposition  of  a  3  x  3 
matrix.  It  is  also  shown  that  generally  there  are  two  distinct  solutions 
for  the  motion  parameters.  Two  results  are  given  in  Part  III.  First,  four 
point  correspondences  between  two  image  frames  are  necessary  and  sufficient 
to  determine  the  eight  pure  parameters.  Second,  with  three  Image  frames, 
the  motion  parameters  are  unique. 
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ABSTRACT 

We  present  a  new  direct  method  of  estimating  the  three-dimensional 
motion  parameters  of  a  rigid  planar  patch  from  two  time -sequential  perspec¬ 
tive  views  (image  frames).  First,  a  set  of  8  pure  parameters  are  defined. 
These  parameters  can  be  determined  uniquely  from  the  two  given  image  frames 
by  solving  a  set  of  linear  equations.  Then,  the  actual  motion  parameters 
are  determined  from  these  pure  parameters  by  a  method  which  requires  the 
solution  of  a  6th-order  polynomial  of  one  variable  only,  and  there  exists  e 
certain  efficient  algorithm  for  solving  a  6th-order  polynomial.  Aside  from 
e  scaling  factor  for  the  translation  parameters  the  number  of  real  solutions 
never  exceeds  two.  In  the  special  case  of  three-dimensional  translation, 
the  motion  parameters  can  be  expressed  directly  as  some  simple  functions 
of  the  8  pure  parameters.  Thus  only  a  few  arithmetic  operations  are  needed. 


I.  INTRODUCTION 


In  the  past,  aost  work  on  motion  estimation  has  been  restricted  to 
two-dimensional  translation.  Recently,  Roach  and  Aggarwal  [1]  and  Huang 
and  Tsai  [2,3,4]  presented  methods  o£  estimating  three-dimensional  motion 
parameters  of  rigid  bodies  based  on  Image-space  shifts.  The  method  of 
Roach  and  Aggarwal  requires  the  solution  of  a  set  of  18  simultaneous  non¬ 
linear  equations;  that  of  Huang  and  Tsai  5  simultaneous  nonlinear  equations. 
Huang  and  Tsai  [4,5]  also  described  a  direct  method  of  estimating  three- 
dimensional  motion  parameters  of  rigid  planar  patches  based  on  the  relation¬ 
ship  between  temporal  and  spatial  differentials  of  image  intensity.  This 
method  results  in  the  solution  of  8  simultaneous  nonlinear  equations.  In 
none  of  the  above  works  was  the  question  of  the  uniqueness  of  the  solution 
to  the  nonlinear  equations  investigated. 

In  this  paper,  we  present  a  new  direct  method  of  estimating  three- 
dimensional  motion  parameters  of  rigid  planar  patches.  We  define  a  set  of 
eight  pure  parameters  and  demonstrate  using  the  theory  of  Lie  Transformation 
Group  that  given  two  pictures,  these  pure  parameters  are  unique.  As  for 
the  estimation  procedure,  we  first  show  using  the  converse  of  the  2nd 
Lie  theorem  [10-13]  that  these  8  pure  parameters  can  serve  as  the  coordinate 
system  of  a  certain  Lie  Transformation  Group.  Then,  we  use  the  result  in 
[10-15]  to  show  that  these  8  pure  parameters  roust  satisfy  a  set  of  linear 
equations.  Furthermore,  the  real  motion  parameters  can  be  computed  from 
these  pure  parameters  by  solving  a  six-order  polynomial. 

Our  new  direct  method  has  several  advantages.  First,  it  requires  the 
solution  of  a  single  sixth-order  polynomial  of  one  variable  only.  Second, 
it  demonstrates  that  more  than  one  solution  may  exist  and  therefore  answers 


the  uniqueness  question.  Third,  in  the  special  case  of  three-dimensional 
translation,  the  motion  parameters  can  be  expressed  directly  as  sane  simple 
functions  of  the  eight  pure  parameters.  Therefore,  only  a  few  arithmetic 
operations  are  needed. 


II.  THE  BASIC  MOTION  EQUATIONS 

We  are  interested  in  estimating  three-dimensional  motion  parameters 
of  rigid  and  deformable  bodies  from  time-sequential  perspective  views 
(frames).  Throughout  this  paper,  we  shall  assume  that  we  work  with  only 
two  frames  at  times  t^  and  t ^  (t^  <  t2>. 

The  basic  geometry  of  the  problem  is  sketched  in  Fig.  1.  Consider  a 
particular  point  F  on  an  object.  Let 

(x,y,z)  -  object-space  coordinates  of  a  point  P  at  time  t.. 
(x',y',z')  *  object-space  coordinates  of  P  at  time 
(X,Y)  -  image-space  coordinates  of  P  at  t^. 

(X',Y')  ■  image-space  coordinates  of  P  at 
It  is  obvious  from  Fig.  1  that 
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Assume  that  from  time  t^  to  t^  the  three-dimensional  object  has  under* 
gone  translation,  rotation,  and  linear  deformation  [8].  Then,  we  have 
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nl  +  n2  +  n3  “  1 

Note  that  (Ax,Ay,Az)  is  Che  amount  of  translation  S  in  the  linear  deforma* 
tion  matrix,  and  (R+I),  where  I  is  a  3x3  unit  matrix,  is  the  rotation 
matrix.  The  rotation  is  around  an  axis  through  the  origin  and  with 
directional  cosines  (n^n^.n^).  The  amount  of  rotation  is  6.  Therefore, 
the  cp^,  cpg*  ?3  defined  in  (5)  are  the  x,  y  and  z  components  of  the  rotation 
vector  with  length  8  and  directional  cosines  (n^.^.n^) . 

Clearly,  Eq.  (2)  represents  an  afine  transformation 
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Conversely,  any  afine  transformation  can  be  decomposed  as  in  Eq.  (2). 


MOTION  OF  PLANAR  PATCHES 


We  now  restrict  ourselves  to  points  on  a  planar  patch  with  equation 

ax  +  by  +  cz  •  1  (8 

at  time  t^.  Then,  it  is  readily  shown  from  Eqs.  (1)  and  (8)  that 

aX  +  bY  +  cF  V 

and  from  Eqs.  (1),  (9),  and  (7)  that 
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We  now  specialize  to  the  case  of  a  rigid  planar  patch.  Then  Eqs.  (11) 
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Eq.  (10)  defines  a  mapping  from  the  2-space  (X,Y)  onto  the  2-space  (X',Y'). 

It  will  be  shown  in  Section  III  that  corresponding  to  any  specified  mapping 
between  the  two  2-spaces,  there  can  be  only  one  set  of  values  for  the 
parameters  a^,...,ag.  We  call  them  the  8  pure  parameters.  In  Section  III, 
we  shall  also  describe  a  method  of  determining  these  pure  parameters  from 
the  two  given  image  frames. 

Once  we  have  determined  the  pure  parameters  a^,a2> . . . ,ag,  we  can 
attempt  to  find  the  actual  motion  parameters:  Ax,Ay,Az .cp^.cp^cp^a.b,  and  c 
by  using  Eqs.  (12).  It  is  obvious  first  of  all  from  looking  at  the  right- 
hand  sides  of  Fqs.  (12)  that  Az  is  a  scale  factor  which  cannot  be  determined. 


We  therefore  let 
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The  unknown  motion  parameters  are  now: 


$>1»  5p2>  ?3»  Ax"»  Ay".  «"»  b",  and  c. 


Thus  we  have  8  nonlinear  equations  with  8  unknowns.  This,  however,  does 
not  mean  that  the  solution  is  necessarily  unique.  In  fact,  it  turns  out 
it  is  not. 

After  some  tedious  manipulations,  we  get  from  Eqs.  (12) 

dgAx"6  +  d5Ax"5  -i-  d^"4  +  d3Ax"3  +  d^"2  -r  djAx"  +  dQ  -  0  (U 


where 
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And  furthermore 
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To  find  Che  motion  parameters,  we  first  solve  Eq.  (14)  for  Ax".  Then 
the  others  are  obtained  from  Eqs.  (17).  Since  Eq.  (14)  is  a  sixth-order 
polynomial  equation,  we  can  have  potentially  6  real  roots  which  give  us  6 
solutions  for  the  motion  parameters.  For  all  the  numerical  examples  we 
have  tried,  only  two  real  roots  are  found  for  Eq.  (14).  One  such  numerical 
exasxple  follows: 
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IV.  DETERMINING  THE  8  PURE  PARAMETERS 

We  now  go  back  and  examine  Eqs.  (10).  For  a  particular  set  of  values 
for  the  parameters  (a^.a^, . . . ,ag) ,  the  equations  represent  a  transformation 
which  maps  the  2-space  (X,Y)  (the  coordinate  space  of  our  image  frame  at 
time  t^)  onto  the  2-space  (X',Y')  (the  coordinate  space  of  our  image  frame 
at  time  t^).  Let  us  consider  the  collection  G  of  transformation  correspond- 

g 

ing  to  all  (a^.a^, . . . ,ag)  6  R  .  We  shall  show  that  it  is  a  continuous  (Lie) 
group  of  dimension  eight  and  that  to  any  given  mapping  from  the  (X,Y)- 
space  onto  (X' ,Y' )-space  corresponds  only  one  set  of  values  for 
(a^ja^, . . . >ag).  Furthermore,  we  shall  describe  a  method  of  determining 
the  pure  parameters  (a^.a^ , . . . ,ag)  from  a  given  pair  of  image  frames  at 
times  t^  and  t 

In  classical  continuous  group  theory,  it  is  known  [13]  that  G  satisfies 
the  four  group  axioms,  namely,  closure,  existence  of  inverse  and  identity, 
and  associativity.  Furthermore,  the  composition  function  for  the  group 
parameters  a^'s  are  continuous.  It  is  also  known  [13]  that  the  a^'s  In  (10) 
are  essential  parameters  in  the  sense  that  the  a^'s  are  functionally  in¬ 
dependent.  However,  it  is  not  known  whether  the  a^s  in  (10)  are  unique,  i.e., 
whether  there  can  be  two  different  sets  of  values  of  a^'s  such  that  (10)  gives 
the  same  mapping  (X,Y)~*  (X'» Y').  Because  of  this  reason,  it  Is  not  easy  to 
verify  whether  G  is  a  Lie  group  according  to  the  modern  definition  since  in 
modern  definition,  in  addition  to  the  properties  satisfied  by  the  classical 
continuous  group  according  to  the  classical  definition  [12,13],  several 
topological  properties  have  to  be  satisfied,  and  these  properties  can  not  be 
easily  verified  unless  we  are  certain  that  the  group  parameters  a^'s  are 
unique.  In  the  following,  we  prove  that  G  is  strictly  a  Lie  group  and  that 
the  a^'s  are  unique. 


Before  we  give  the  proof,  we  motivate  it  by  the  following  considera¬ 
tions.  Let  us  assume  that  G  is  indeed  a  Lie  group  and  that  (a^^, . .  •  ,ag) 
is  a  coordinate  representation  for  the  group  G.  The  identity  element  of 
the  group  is  obviously  e  *  (1,0, 0,0, 1,0, 0,0) .  Then  the  operators  of  the 
Lie  algebra  associated  with  the  group  G  are  given  by: 
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where  g  is  used  to  represent  a  member  of  G.  From  Eqs.  (10)  we  get  readily 
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Now  we  stare  our  proof.  Consider  the  set  of  vector  fields  on  the 
differentiable  manifold  (X,Y)  as  given  by  Eqs.  (20).  It  can  be  easily 
verified  that  none  of  the  X,  can  be  expressed  as  a  linear  combination  of 
the  others,  i.e.,  J-i,2,...,8}  are  linearly  independent;  and  furthe 

more  for  any  1,* 
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where  ere  constants.  From  these  two  properties  of  {X^}  we  conclude 
from  the  converse  of  Lie's  second  fundamental  theorem  [10-13]  that  there  is  a 
unique  Lie  group  of  transformation  of  order  8  which  has  fx^}  as  its  Lie 
algebra.  We  proceed  now  to  show  that  G  is  that  group. 

From  [15,10,11] ,  one  can  generate  the  finite  equations  of  a  Lie  group 
of  transformation  with  the  X^s  as  the  canonical  coordinates  of  the 
second  kind  as  follows: 


X'  •  exp&gXg)exp ...  exp&jXpx 

Y'  *  exp (X gXg )exp (kjXj)  ...  expfl.jX^Y 

It  is  proved  in  [6]  that  (22)  is  equivalent  to  the  following: 
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Comparing  (10)  and  (23)  shows  that  G  is  indeed  a  Lie  Group  of  trans¬ 
formation,  and  that  since  the  X^'s  are  the  canonical  coordinates,  they  are 
unique,  and  therefore  from  (24),  the  pure  parameters  a^  ...agare  also  unique. 


Let  f  be  any  function  defined  on  the  2-space  (X,Y)  (in  our  case,  f  will 
be  the  intensity  of  the  picture  elements),  then  from  Lie  Group  Theory 
[10-15],  we  have 
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Af  -  Z  3.X  f  (25) 

i-1  1  1 

where  ^i  "  ai  “  ei 

e^  »  ith  component  of  the  group  parameters  at  the  identity 

Af  -  f(X*  ,Y')  -  f(X,Y)  ■  frame  difference 

(The  implicit  assumption  here  is  that  the  intensities  of  the  picture 
elements  in  the  image  frames  corresponding  to  the  same  physical  object 
point  are  the  same.)  Clearly, 

*!«&!+  1 
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®t  •  01,  i  ■  2 ,3,4,6 ,7 ,8 

Eq.  (25)  is  used  to  determine  the  g^’s  and  therefore  the  8  pure  parameters 
a^'s.  We  pick  8  or  more  points  (X,Y),  calculate  at  each  point  Af  and 
Xjf  (i  ■  1,2,..., 8),  and  substitute  into  Eq.  (25)  to  obtain  8  or  more 
equations  which  are  linear  in  the  8  unknowns  g^’s.  Then  we  find  the  least- 
square  solution. 

V.  DISCUSSIONS 

la  this  paper  vs  have  investigated  ehe  problem  of  estimating  three- 
dimensional  motion  parameters  of  a  rigid  planar  patch  from  two  image  frames 
The  following  results  have  been  established: 
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1)  The  face  chat  we  can  define  8  pure  parameters  a^a^...^^  which  ara 

unique  foe  any  given  napping  from  eha  2 -space  (X,?)  (the  coordinate 
space  of  the  inage  frame  at  t1 )  onto  the  2-space  (X',Y')  (the  coordinate 

space  of  the  image  frame  at  c^). 

2)  A  method  of  determining  Che  actual  moeion  parameters  p^,  p^,  "» 

iy",  a",  b",  and  c  from  eha  pure  parameters  a^,  a^,  ....  a^  which 

requires  the  solution  of  a  6th-order  polynomial  of  one  variable 

only.  Aside  from  a  scaling  factor  in  the  translation  parameters, 
the  number  of  real  solutions  never  exceeds  two. 

3)  A  method  of  determining  the  3  pure  parameters  a^,  a^ ,  . ..,  a^  from 
the  two  given  image  frames.  This  requires  the  solution  of  a  set  of 
linear  equations  only. 

It  is  to  be  noted  chat  1)  and  2)  are  Independent  of  3).  The  pure 
parameters  can  be  determined  by  ocher  methods.  For  example,  if  one  can 
identify  8  or  more  corresponding  point  pairs  in  Che  two  image  frames  [2]» 

Chen  the  a^'s  can  be  determined  from  Eqs.  (10)  by  solving  a  set  of  linear 
equations. 

Recently,  an  alternative  way  of  analyzing  the  uniqueness  problem  and 
estimating  the  three-dimensional  motion  parameters  has  been  developed  which 
stems  from  the  results  contained  in  this  paper,  and  requires  computing  the 
singular  value  decomposition  (SVD)  of  a  certain  3x3  matrix  only.  The  eight 
pure  parameters  defined  by  the  authors  in  this  work  will  again  be  used. 

It  is  briefly  mentioned  in  [6 ] ,  and  the  detailed  paper  will  be  submitted 


soon. 


Furthermore,  the  conclusion  that  the  motion  parameters  are  generally 


not  unique  is  of  course  independent  of  the  method  of  determining  these 
parameters.  The  question  arises:  Are  the  motion  parameters  unique,  aside 
from  the  scaling  factor,  if  the  rigid  patch  is  nonplanar?  We  have  solved 
this  problem  recently  [7],  and  will  publish  it  in  the  near  future. 
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Abstract 

Wa  show  that  the  three-dimensional  motion  parameters  of  a  rigid  planar 
patch  can  be  determined  by  computing  Che  singular  value  decomposition  (SVD) 
of  a  3  x  3  matrix  containing  the  eight  so  called  "pure  parameters".  Further¬ 
more,  aside  from  a  scale  factor  for  the  translation  parameters,  the  number  of 
solutions  is  either  one  or  two,  depending  on  the  multiplicity  of  the  singular 
values  of  the  matrix. 


I.  INTRODUCTION 


The  processing  of  image  sequences  involving  motion  has  become  increas¬ 
ingly  important.  Because  of  the  key  role  motion  estimation  plays  in  image 
sequence  processing,  a  considerable  amount  of  effort  has  been  devoted  to  this 
topic,  for  example,  see  Ref.  [1-17].  However,  except  for  [1-8] [16] [17] ,  most 
past  work  considers  only 2 -D  motion,  especially  translation.  Ref.  [1-8] [16] [17] 
were  among  the  first  to  consider  3-D  motion  and  [1][3][17]  were  among  the 
first  to  consider  .the  problem  of  uniqueness  of  solutions.  In  [1][3],  the 
eight  "pure  parameters"  were  introduced  for  the  case  of  a  rigid  planar  patch 
undergoing  general  3-D  motion,  and  proved  to  be  unique  given  two  successive 
(in  time)  perspective  views.  The  proof  makes  use  of  the  theory  of  Lie  Group 
of  transformations.  It  .was  also  shown  that  these  8  pure  parameters  can  be 
computed  by  solving  a  set  of  linear  equations.  Furthermore,  once  the  pure 
parameters  are  determined,  the  actual  motion  parameters  can  be  computed  by 
solving  a  6-th  order  polynomial  equation  of  one  variable  if  the  motion  is 
small.  Theoretically,  the  number  of  solutions  cannot  exceed  six  aside  from 
a  common  scale  factor  for  the  translation  parameters;  experimentally,  the 

number  of  solutions  has  been  found  to  be  two.  In  this  paper,  we  show 
chat  whether  the' motion  is  small  or  not,  once  the  eight  pure  parameters  are 
computed  using  the  method  described  in  [1][  3],  Che  actual  motion  parameters 
can  be  estimated  by  computing  the  SVD  of  a  3  x  3  matrix  consisting  of  the 
eight  pure  parameters.  Also,  by  using  the  rigidity  constraint  and  the  fact 
that  a  plane  in  3-space  can  be  oriented  in  at  most  two  possible  ways  in  order 
to  intercept  an  ellipsoid  at  a  circular  cross-section,  we  prove  that  the 
number  of  solutions  is  either  one  or  two,  depending  on  the  multiplicity  of  the 
singular  values.  Physical  description  of  the  motion  is  stated  and  justified. 


II.  THE  EIGHT  FUSE  PARAMETERS  AND  THE  MOTION  OF  PLANAR  PATCHES 
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The  basic  geomecry  of  Che  problem  considered  in  [I] [3]  is  repeated 


here  in  Fig.  1.  Throughout  this  paper,  we  shall  assume  that  we  work  with 
only  two  frames  at  time  t^  and  t^  (t^  <  t2).  Consider  a  particular  point  P 
on  the  object.  Let 


<x,y,z)  - 
(x'.y'.z*)  - 
(X,Y)  - 
(X ' ,  Y ' )  - 


object-space  coordinates  of  a  point  P  at  time  t^. 
object-space  coordinates  of  P  at  time  t2< 
image-space  coordinates  of  P  at  t.. 
image-space  coordinates  of  P  at  t^. 


It  was  shown  in  £  1] [31  that  for  a  rigid  planar  patch  undergoing  3-D 
motion  (a  rotation  with  a  small  angle  0  around  an  axis  through  the  origin 
with  directional  consinas  n^  n2>  n^  followed  by  a  translation  with  trans¬ 
lation  vector  (Ax,  Ay,  Az)) .the  image-space  coordinates  before  and  after  the 
motion  are  related  by  the  following  equations: 


«jX  +  a2Y  +  a3 

•t*  +  V+  1 


a4X  +  «3Y  +  «6 
»7*  +  »8Y  +  1 


where 


1  +  a-Ax 
1  +  c«Az 

-n^S  +  b*Ax 

1  +  c-Az 

(n20  +  c*Ax)F 

1  +  c*Az 

n30  +  a* Ay 

“l  +  c*Az 


l  +  b*Ay 
1  +•  c»Az 
(-n^0  +  c.Ay)F 
l  +•  c*Az 
-n20  +  a*Az 
(1  +  c*Az)F 
n^0  +  b«Az 
(1  +  c*Az)F 


I 


Fig*  1  Basic  g soma cry  fot  ehraa- 
diaaasionsL  noeiott 
as ciaac ion. 


where  a,  b,  e  are  the  parameters  chat  appear  la  the  following  equation: 


ax  +  by  +  cz  -  1 


which  describes  the  surface  of  the  object  in  the  object  space  coordinate 
system  at  time  t^. 

Eqs.  (2)  are  applicable  when  Che  rotation  angle  is  small.  For  general 
3-0  motion,  it  is  well  known  in  mechanics  and  computer  graphics  [18]  that  any 
3-0  rigid  body  motion  is  equivalent  to  a  rotation  followed  by  a  translation 


R  y  Ay 


where  R  is  a  3  x  3  orthonormal  matrix 


?  2 
n^+Cl-n^cosG 

nln2  ^  *"cos® 

n^n  ^  ( l-cos9  )  -n2  s  in9 


njttj (l-cos0)-njSin9  njn^l-cosOJ+njSinO 


2  2 
n|+(l-n|)cos9 


n2n2(l-cos0)-nj8ia9  (5 


2  2 

njn^Cl-cosgj+n^sing  n^+Cl-npcosO 


Following  exactly  the  same  procedure  as  in  [1][3],  one  can  show  that 


(1)  is  again  valid  if  (2)  is  replaced  by 
2  2 

nt  +  (l-nr)cos8  +  a. Ax 


n.n. (l-cos0)  —  n,sln8  +  b-Ax 


nr  +  (l-n£)cos0  +  c*Az 


+  (l-n2,)cos0  +  b*Ay 


n|  +  (l-n^)cos0  +  c*Az 
n,n-(l-cos0)-n„sin0+a* Az 


+  (l-n£)cos0  +  c*Az 
where  for  simplicity  we  have  set  F 


n^  +  (l-n|)cos0  +  c*Az 
[n-n-(l-cosS)+n.  sin0+c*Ay]F 


n^  +  (l-n|)cos0  +  c*Az 
n-n-Cl-cos©)^,  sin0+b*  Az 


n%  +  (l-n4)cos0  +  c*Az 


K 
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It  was  shown  in  [1JI3J  using  Lie  Group  Theory  that  given  two  perspective 
views  at  t ^  and  t2,  the  eight  pure  parameters  a^,  a2,  ...  ,  ag  are  unique, 
and  they  can  be  estimated  by  solving  a  system  of  linear  equations. 


2- 


III.  COMPUTING  ACTUAL  MOTION  PARAMETERS  FROM  PURE  PARAMETERS 


By  using  the  rigidity  constraint  and  the  fact  that  a  plane  in  3-space 
can  be  oriented  in  at  most  two  possible  directions  in  order  to  cut  a  circle 
in  an  ellipsoid,  we  shall  prove  that  the  number  of  possible  solutions  for 
the  motion  parameters  can  never  exceed  two  aside  from  a  scale  factor  for  the 
translation  parameters.  The  number  of  solutions  depends  upon  the  multiplicity 
of  the  singular  values  of  the  following  matrix  A  consisting  of  the  eight 


pure  parameters  a^'s : 


a4  a5 


The  SVD  of  A  is  given  by 


V  -  UAV 


where  X^'s  are  the  singular  values  of  A,  and  U,V  are  3x3  orthonormal 


matrices. 


Let  k  ■  n2  +  (l-n?)cos9**  cdz;  then  it  can  be  readily  shown  that 


A  "  k 


kA  - 


I# 


From  (3)  and  (4),  it  can  be  seen  that 


-  _ 

»  * 

r  ,  i 

•  « 

x’ 

X 

Ax 

£a  b  cj 

X 

y’ 

-  R 

y 

+ 

Ay 

y 

z * 

»  - 

z 

■  m 

Az 

P  0 

z 

.  * 

/ 

►  < 

*  « 

Ax 

t* 

b 

C1 

X 

W 

R  + 

Ay 

y 

Az 

z 

•  « 

4 

»  « 

-  k  A 


x 

y 

z 


(10) 


If  we  transform  the  original  coordinate  system  with  the  orthonormal  matrix 


V  in  (8)  as  depicted  in  Fig.  2,  where  (xq,  yQ,  zq)  is  the  new  coordinate 


system  after  transformation,  we  have 


and 


(U) 


(12) 


Substituting  (11),  (12)  into  (10)  gives 
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yn 

Fig.  2  The  relationship  between  the  (x,y,z)  and 
(xn,yn,zQ)  coordinate  systems. 


•  m 

■  e 

x' 

X 

n 

n 

K 

-  k  A  V 

yn 

z' 

z 

n 

■  ■ 

n  . 

By  taking  the  Euclidean  norms  of  the  vectors  on  both  sides  of  (13) 
we  obtain 


H  yn  Vl  V 

»  e 

x' 

n 

*  k2  [\  ■.]  vTaI*v 

•  « 

X 

n 

K 

•i 

S  < 

.v 

T 

Since  V  is  orthonormal,  V  *V  on  the  left  hand  side  can  be  replaced 
by  an  identity  matrix.  Substituting  (8)  into  (14)  gives 


Replacing  U  U  In  (IS)  by  an  identity  matrix  gives 


k2(Xlxn  +  AA  + 


(16) 


This  is  the  key  equation  that  will  lead  us  to  the  solution  o£  the 
uniqueness  problem,  as  will  be  seen  hereafter. 

In  the  following,  we  state  and  prove  three  theorems  regarding  the 
uniqueness  and  computation  of  the  motion  parameters  given  the  pure  parameters, 
and  the  physical  characterizations  of  the  motion  in  the  object  space  for 
different  multiplicities  of  the  singular  values  of  the  matrix  A. 


THEOREM  I 


If  the  multiplicity  of  the  singular  values  of  A  is  two,  e.g.,  X.,  ■ 

X2  t  X3,  then  the  solution  for  the  motion  and  geometrical  parameters  is  unique 
aside  from  a  common  scale  factor  for  the  translation  parameters,  and 


-1  a3  T 

R  -  XL  A  -  -  s)  U3  V3  , 


Ax 

Ay 

Az 


-1  X3 

»  l(^  -  .)  V3 


and 


w  V, 


where  s  »  det(U)det(V) 

w  is  a  scale  factor 

a,b,c  are  the  parameters  in  (3)  which  is  the  planar  equation  of  the 


object  surface  at  time  t^ 

U„  ,V,  are  the  third  columns  of  U  and  V  respectively. 

Furthermore,  a  necessary  and  sufficient  condition  for  the  multiplicity  of 
the  singular  values  of  A  to  be  2  is  that  the  motion  can  be  realized  by 
rotating  the  object  around  the  origin  and  then  translating  it  along  the  normal 
direction  of  the  object  surface. 


[Proof] 

The  two  sides  of  (16)  can  be  equated  to  a  collection  of  positive 
values  corresponding  to  the  range  of  values  for  x,  y,  z  and  x',  y',  z'.  Let 
d  be  one  such  value.  Then  we  have 


and 


.,2 


+  y;2  +  <2 


.2„2  2  ,  .2  2  .2  2. 

-  k  (\Lxa  +  X2yn  +  \2za) 


(17) 

(18) 


Clearly,  (17)  defines  a  sphere  in  the  (x^,  y^,  z^)  space,  while  (18) 
defines  an  ellipsoid  la  the  (xq,  yQ,  zQ)  space.  Since  X^  -  X^,  two  of  the 
three  principal  axes  of  the  ellipsoid  are  equally  long.  Since  the  object 
surface  is  assumed  to  be  planar,  the  collection  of  the  points  on  the 
object  surface  that  also  satisfy  (17)  must  be  the  circle  which  lies  on  the 
Intersection  of  the  sphere  and  the  object  surface  at  time  t^  (see  Fig.  3). 
Because  of  (16),  (17)  and  (18),  all  the  points  on  this  circle  at  time  tj 
must  also  satisfy  (18)  at  time  t^,  l.e.  they  oust  lie  on  the  Intersection 
of  the  object  surface  and  the  ellipsoid.  Due  to  the  rigidity  constraint. 


this  intersection  should  also  be  a  circle.  But  the  only  possibility  for  a 
plane  to  cut  a  circle  out  of  an  ellipsoid  with  two  of  the  three  principal 


(d) 
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Cb) 

Fiq.  3  (a)  .  The  ODject  surface  always  intercepts  the  sphere  defined 

by  (17)  at  a  circular  cross-section. 

(b)  .  Tne  ooject  surface  Intercepts  the  ellipsoid  defined  by 

(ltf)  at  a  circular  cross-section  only  if  it  is  properly 


oriented 


axes  equally  long  Is  that  Che  plane  be  perpendicular  Co  Che  major  axis  (Che 
longesC  one)  of  Che  ellipsoid,  as  depicCed  in  Fig.  4.  Thac  is  Co  say,  Chere 
is  only  one  possible  orlencacion  for  Che  objecC  surface  before  Che  mocion. 
This  is  Che  key  sCep  chat  will  lead  us  Co  Che  conclusion  chac  Che  mocion 
parameters  are  unique,  as  will  be  seen  shortly. 

Note  chat  X^  (and  X^)  caa  never  be  zero  since  were  Chis  Co  be  true, 

Che  ellipsoid  defined  by  (18)  would  have  degenerated  into  Cwo  parallel  planes 
and  chere  is  no  way  Che  objecC  surface  can  intercept  cwo  parallel  planes  at 
a  circular  cross-section. 

Zn 

y 
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Fig.  4.  The  object  surface  muse  be  perpendicular  Co  Che  z  axis 
Co  intercept  the  ellipsoid  at  a  circular  cross-section 
when  X.  ■  X-  +  X,. 


Since,  as  depicted  in  Fig.  4,  the  object  surface  must  be  perpendicular 


to  the  z  axis,  and  since  the  z  axis  is  obtained  by  rotating  the  z  axis 
n  n 

with  the  orthonormal  matrix  V  as  in  Fig.  2,  it  is  seen  that 


a 

b 


0 

0 


w  V, 


(19) 


where  a,b,c  are  the  parameters  in  (3), 

Vj  is  the  third  column  of  V  in  (8),  and 
w  is  an  arbitrary  constant. 


Substituting  (19)  and  (8)  into  (9)  gives 


R  + 


r  i 

Ax 

Ay 

Az 


[o  0  w]  •  v’ 


k  UAV 


(20) 


Premultiply  (20)  by  U  and  postmultiply  (20)  by  V  to  give 


T  T 

0  R  V  +  0 


Ax 

Ay 

Az 


1 0  0  w 


]  VTV 


kA 


or 


_ 

r  n 

R'  + 

Ax' 

Ayi 

Az' 

•  m 

M 

■ 

T* 

o 

o 

X1 

X2 

X3 

(21) 


(22) 


where 


R'  -  UTR  V 


(23) 


2 


*  « 

Ax' 

Ax 

Ay' 

T 

»  U 

Ay 

Az' 

» 

L  Lz  . 

(22)  gives 


R' 


kXj^  0  -w*Ax' 
0  kX^  -w*Ay' 
0  0  IcXj-w-Az1 


(24) 


(25) 


It  will  be  shown  now  that  Ax'  and  Ay'  in  (25)  are  zero,  therefore  R' 
is  diagonal. 

Since  U,  V  and  R  in  (23)  are  all  orthonormal,  R'  is  also.  Taking  the 
inner  product  o£  the  2nd  and  3rd  columns  of  R',  and  equating  it  to  zero  gives 


kXjW  •  Ay'  ■  0  (26) 

X^  and  k  cannot  be  zero  since  were  X^  or  k  to  be  zero,  the  1st  and  2nd 
columns  of  (25)  would  be  zero,  which  contradicts  the  fact  that  R'  is 
orthonormal.  Obviously,  w  cannot  vanish  either,  otherwise,  a,  b  and  c  would 
vanish,  which  contradicts  (3).  Therefore,  (26)  implies  that  Ay'  ■  0. 

Similarly,  one  can  show  that  Ax'  ■  0.  Thus,  (25)  is  diagonal,  which, 
when  combined  with  the  fact  chat  R1  is  orthonormal,  gives  the  following: 


kX^  •  +1  or  -1 

kXj  -  wAz'  »  +1  or  -1 

Ha  show  that  k  has  to  be  positive: 

From  (10),  we  have 

*'  •  fc(a?x  +  V  +  *> 


(27) 

(28) 


(29) 
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t-. 


For  x  -  0  and  y  ■  0,  z*  ■  kz.  Since  the  object:  muse  be  in  front  of 
Che  camera,  z  and  z'  are  both  positive,  which  implies  that  k  is  positive. 

Since  is  nonnegative  by  definition,  the  right  hand  side  of  (27)  can 
not  be  -l.  Therefore, 


and 


(30) 


where  s  is  either  +L  or  -l. 

Since  R'  ■  URVT,  we  have  det(R')  ■  det(U)det(R)det(V) 
Thus  a  -  det(U)det(V).  (28)  gives  the  following: 

dz'  ■  w”^(r^  -  ») 

From  (24),  (31)  and  the  fact  that  dx'  *  dy'  ■  0,  we  have 


det(U)det(V). 


(31) 


k  V 


►  « 

dx 

*  m 

dx' 

0 

dy 

-  U 

dy* 

-  U  • 

0 

dz 

dz' 

-1  *3 

w  A(r®  -  a) 

•  « 

L  \ 

-l  *3 

w  l(^  -  s)  U3 


(32) 


(19),  (20)  and  (32)  imply  that 


R  -  X^A  - 


dx 

dy 

dz 


2 


-1  3  T 

H  A  -  -  *>  "j  v3 


(33) 


In  the  following,  it  will  be  shown  that  the  solutions  for  the  rotation 
matrix  R  in  (33)  is  unique,  and  that  aside  from  a  scaling  factor,  the  trana- 
lation  parameters  Ax,  Ay,  Az  in  (32)  and  the  geometrical  parameters  a,  b,  c 
in  (19)  are  also  unique. 

The  first  thing  to  show  is  that,  once  A  is  given,  is  fixed  except 
for  the  sign.  From  (8), 


T  2 

A  A  U3  -  \2  U3 


(34) 


Let  Q  be  any  orthonormal  eigenvector  matrix  of  A*A  .  Then 


A  A1  -  QT 


(35) 


From  (34)  and 

(35) 

,  we  have 

/ 

« 

1 _ 

qt 

4 

.  2 

Q  '  X3  1  ? 

* 

X3 

4 

or 


PDj  >  0 


where 


A  T  2  1 

P  -  AA  -  X3I  -  Q 


2  2 
X  -X 
*1  j 


2  2 
X  -X* 
*1  *3 


(36) 


P  has  rank  2  since  X^  -  X3  on  the  diagonal  of  the  diagonal  matrix  in 


T  2 

(36)  is  nonzero.  Also,  F  is  fixed  once  A  is  given  since  P  *  AA  -  X‘l 
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Therefore,  U3  is  fixed  except  for  the  sign. 

The  next  thing  to  prove  is  that  (33)  is  unique,  i.e. 

Ax  £a  b  cj 
Ay 
Az 

is  unique  once  A  is  given.  Two  cases  are  to  be  discussed.  The  first  is 

when  X3  1*  0.  In  this  case,  s  -  sgn(det(A))  since  A  •  U  A  V  and  thus 

s  »  det(U)det(V)  ■  (XjX2X3)’*det(A)  ■  sgn(det(A)).  Therefore,  given  A, 

T 

s  is  fixed.  The  next  thing  to  prove  is  that  in  (33)  is  unique. 

Since  U3  and  are  fixed  except  for  the  sign,  all  one  has  to  show 

is  that  when  changes  its  sign,  oust  also. 

From  (8)  we  have 


A  V  -  U  A  -  [Xft  >.2U2  A3D3] 

thus 

A  V3  -  X3U3 

Since  A  and  X3  are  fixed  given  two  perspective  views,  we  see  that  when 

T 

V3  changes  its  sign,  U3  must  also.  Therefore,  U3V3  has  fixed  sign.  We  have 
thus  proved  that  the  product 


and  therefore  R  in  (33)  are  all  unique. 

For  Che  second  case  for  which  X3  ■  0,  we  have  from  (33)  that 
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R  -  X'lA  +  S  U3V3 

-  X’LA  +  det(U)  •  U3  •  det(V)  •  V*  (37) 

l£  U3  changes  its  sign,  det(U)  will  also.  Thus  the  sign  of  det(U)  in 

T 

(37)  remains  unchanged.  Similarly,  the  sign  of  det(V)  V3  also  is  fixed  when 

V3  changes  its  sign.  Therefore,  the  uniqueness  of  (37)  is  not  shaken  by 

the  ambiguity  of  the  signs  of  U3  sad  Vy 
T 

Since  A  A  have  double  eigenvalues,  the  eigenvectors  and  that 
correspond  to  the  multiple  eigenvalues  X^  (  ■  X^)  are  orthonormal  to  each 
other  but  may  be  anywhere  in  a  certain  fixed  plane  perpendicular  to  V3  (Note 
that  we  are  now  interpreting  eigenvectors  geometrically  as  some  vectors  in 
3-space.)  If  the  order  of  and  V2  on  the  plane  are  interchanged  while 
keeping  V3  fixed,  the  sign  of  det(V)  will  change.  We  are  now  to  prove  that 
when  this  does  happen,  the  sign  of  det(U)  will  also  change,  thereby  keeping 
(37)  fixed.  It  is  obvious  from  (8)  that 

(X^A)  VL  -  UL 

(X’lA)  V2  -  U2 

Since  \1  and  A  are  fixed,  when  and  V2  are  interchanged,  and  U2 
will  also.  Therefore,  when  det(V)  changes  sign,  det(U)  will  also.  Thus  for 
the  case  when  X3  “  0,  Che  product 


as  well  as  R  in  (33)  are  unique . 


We  are  now  to  prove  that  a  necessary  condition  for  ■  X2  J1  X^  is 

that  the  translation  vector  Ax  is  parallel  to  the  normal  direction 

Ay 


of  the  object  surface  at  time  t2- 

Since  before  the  motion,  i.e.  at  t^,  f"  a  1  is  normal  to  the  object 


surface,  this  vector  is  rotated  by  R  at  time  tg  and  becomes  R  b 


It  is  only  necessary  to  prove  that  there  exists  a  scalar  q  such  that 


q  R  b 


From  (19),  (32),  (33)  and  (38),  we  have 

-  s>  °3  -  •  a>  u3  “  V3 


<1  w3.'lA  v3  -  ,  „  .  ,)  #3  V*  v3 


A  *  Xj  +  Aj  ^2  Vg  +  \3  Uj  V*  f  thus  A  V3  in  (39)  becomes 


A  V3  "  **1U1  V1  V3  +  ^2U2  V2  V3  +  *3U3  V3  V3 


0  + 


0  +  X3U3 


Substituting  (40)  into  (39)  gives 


-1  ^3  ^3  ^3 

w  A(r-  -  •)  »,  -  q  w  (r—  -  ~  +  s)  u, 


q  w  s  U, 


Therefore,  if  we  take  q  as  w  s  (? -  s) ,  then  (38)  will  be  satisfied.  We 

A1 

have  thus  proved  that  the  necessary  condition  for  X^  »  X2  +  X^  is  that  the 
motion  can  be  realized  by  first  rotating  the  object  around  an  axis  passing 
through  the  origin,  and  then  translating  it  along  the  normal  direction  of 
the  object  surface  at  time  t^> 

In  the  proof  of  Theorem  II,  it  will  be  shown  that  if  the  translation 
is  along  the  normal  direction  of  the  object  surface  at  t^>  then  the  singular 
values  of  A  can  not  be  all  distinct.  This  fact,  together  with  Theorem  III, 
provide  the  sufficiency  part.  Q.E.D. 


THEOREM  II 


If  the  singular  values  of  A  are  all  distinct,  e.g. ,  X1  >  X^  >  *3* 
then  there  are  exactly  two  solutions  for  the  motion  and  geometrical  param¬ 
eters  aside  from  a  scale  factor  for  the  translation  and  geometrical  param¬ 
eters,  and  that 

0 
1 
0 


8 

0 

sa 


Ax 

Ay 

Az 


X3 

(rr  -  sa)  uj 

A2  ^ 


a 

b 


w  (6V.  +  V,] 


where  « 


xl  -  X* 

2  3 


X1  +  sX36a 
X2(l  +  52) 


8  -  +/1  -  ot* 


s  -  det (U) dec (V) 


(in  each  of  Che  two  solutions,  sgn(B)  ■-sgn(5).) 

Furthermore,  a  necessary  and  sufficient  condition  for  distinct  singular 
values  is  that  the  motion  can  be  decomposed  into  rotation  around  an  axis 
through  the  origin  followed  by  translation  along  a  direction  different  from 
the  normal  direction  of  the  object  surface  at  time  t^. 

[Proof] : 

Since  the  three  singular  values  are  distinct,  the  three  principal 
axes  of  the  ellipsoid  defined  by  (18)  are  of  different  lengths.  By  using 
the  same  argument  as  in  Theorem  I,  the  object  surface  at  t^  must  be  oriented 
in  such  a  way  that  it  cuts  a  circle  out  of  this  ellipsoid.  It  is  easy  to 
verify  using  basic  analytical  geometry  that  a  plane  can  be  oriented  in  only 
two  possible  directions  (see  Fig.  3)  in  order  to  cut  a  circle  out  of  an 
ellipsoid  whose  three  principal  axes  are  of  different  lengths. Since  X1  >  X2  >  X 
,the  longest  principal  axis  is  aligned  with  zn  axis, and  the  vector  normal  to 
the  object  surface  is 


0  ,  where  6  ■  ± 


in  the  (x  ,  y  ,  z  )  coordinate  system, and  VI  Is  some  constant, 
n  n  n 


FP-rzrr 


"u 


1 

Fig.  5  The  object  surface  can  be  oriented  in  exactly  two  possible 
directions  in  order  to  Intercept  an  ellipsoid  defined  by  Cl 
at  a  circular  cross-section  wnen  ^,>  A*  >  A*  • 
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Since,  as  shown  in  Fig.  2,  the  (x  .  y  ,  z  )  axes  are  obtained  by 

n  n  n 

the  (x,  y,  z)  axes  with  the  orthonormal  matrix  V,  we  see  that 


•  - 

r*  - 

a 

£ 

b 

-  w  V 

0 

c 

1 

-  _ 

where  w  is  some  constant.  Substituting  (41)  into  (9)  gives 


kA  *  k  U  A  V  -  R  +  w 


Ax 

Ay 

Az 


[*  0  *1 


1  V 


Premultiplying  (42)  by  UT 


and  postmultiplying  (42)  by  V  give 


kA  -  UT  R  V  +  w  UT 


Ax 

Ay 

Az 


Thus 


[*  0  lJ 

where  R'  4  0T  R  V 


Ax' 

Ax 

Ay’ 

A  w  UT 

Ay 

Az' 

•  « 

Az 

•  m 

R'  -  k  A  - 


Ax' 

Ay' 

Az' 


l]  VT  V 


rotating 


(41) 


(42) 


(43) 

(44) 


(45) 


From  (43) , 
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kX^fi'Ax' 

0 

-Ax' 

-  6*  Ity9 

kX2 

-Ay' 

-5‘Az ' 

0 

kX3~Az 

(46) 


Since  U»  R  and  V  are  orthonormal,  (46)  implies  that  R'  is  also  orehonormal. 
Taking  the  inner  produce  of  columns  2  and  3  and  equating  the  result  to  zero 
gives 

kX2  •  Ay'  -  0  (47) 

Since  >  *2  >  *3  —  0,  we  have  >  0.  Therefore,  (47)  implies  that 
Ay'  ■  0.  Thus  (46)  becomes 


R' 


'kXj-6-Ax'  0 

0  kX2 

-fi-Az'  0 


-Ax' 

0 

kX3-Az* 


(48) 


The  normality  of  column  2  implies  that  kX2  •  ±  1.  But  since  k  >  0  and 
X2  >  0,  we  have  kX2  ■  1,  or  k  ■  l/x2*  Furthermore,  from  the  fact  that 
columns  1  and  3,  as  well  as  rows  1  and  3  of  R'  are  mutually  orthogonal,  and 
that  the  norms  of  the  rows  and  columns  of  R'  are  unity,  it  can  be  shown  that 


where 


^1 

o  ■  - 6- Ax'  -  s(kX3  -  A*') 


g  ■  -Ax' 


sfi«A*’ 


(49) 


(50) 

(51) 


8  ■ 


det(U)det(V) 


(52) 
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Since  U  and  V  are  orthonormal,  from  (52),  s  Is  either  +1  or  -1.  It  is  to 
be  shown  that  although  det(U)  and  det(V)  may  be  +1  or  -1  for  a  particular  A, 
s  is  unique  once  A  is  given. 

T 

Recall  that  U^,  and  ere  the  eigenvectors  of  AA  corresponding  to 
2  2  2 

eigenvalues  X^,  X2»  and  X^  respectively.  Since  X^,  X^  and  X^  are  distinct, 
U1#  U2,  ^3*  V]_»  V2  V3  ere  all  fixed  except  for  the  signs.  However,  as 
was  seen  in  the  proof  of  Theorem  I,  we  have 


A  V 

A  V 

A  V 


1 

2 

3 


X1U1 


X2U2 


x3u3 


Therefore,  when  changes  Its  sign,  will  also,  where  i  -  1,  2,  3. 
Hence  the  sign  of  det(U)det(V)  remains  fixed.  Thus,  s  is  unique. 
From  (50)  and  (51) ,  we  have 


*1 

a  -  -  8-6  (53) 

A2 


*3 

6(a+  s~)  -  6  (54) 

a2 


Cancelling  8  in  (53)  and  (54)  gives 

X1  +  sx362 

a  »  - r 


x2d  +  o 


where 


6  - 


From  (50)  end  (51),  end  the  fact  that  Ay'  ■  0,  we  have 


m  m 

•  - 

Ax' 

-e 

ay* 

m 

0 

X3 

Az' 

—  -  sc r 

«•  « 

L  X2 
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(55) 


Fro®  (45),  (47)  and  (55), 


m  a 

Ax 

Ax' 

-e 

Ay 

■  w"1  U 

Ay’ 

m  W-LU* 

0 

Az 

Az* 

X3 

T"  *  *« 

■  a 

m  m 

-  £ 

w'1!-®  UL  +  (—  -  sat)  U3J 


(56) 


From  (44)  and  (49), 


R  -  U 


a  0.0 
0  l  0 

J.-80  0  8  Of 


(57) 


Fro®  (41),  (56),  (57),  and  Cha  fact  that  s  Is  fixed,  we  see  that  there 
are  exactly  two  solutions  aside  fro®  a  scaling  factor  for  the  translation 
and  geometrical  parameters. 

It  is  to  be  shown  that  a  necessary  and  sufficient  condition  for  the 
singular  values  to  be  distinct  is  that  the  translation  vector  is  not  aligned 
with  the  normal  direction  of  the  object  surface  after  rotation  (or  at  time 
tj).  The  sufficiency  part  was  proved  in  Theorem  I.  The  necessity  part  is 
proved  by  contradiction.  We  shall  show  that  if  the  translation  vector  is 
along  the  normal  direction  of  the  object  surface  at  tj,  then  the  singular 


f 


values  cannot  be  distinct. 

It  was  Indicated  in  the  proof  of  Theorem  I  chat  Che  normal  direction 


c 


r  i 

a 

Ax 

of  the  object 

surface  at  t.  is  aligned  with  R 

b 

.  Suppose 

Ly 

•  m 

| 

c 

A* 

a 

a  * 

L 

parallel  to  R 

b 

,  then 

c 

i 
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is 


Ax" 

a 

Ly 

-  h  R 

b 

Lz 

W  m 

c 

m  .J 

(38) 


for  some  constant  h.  With  (41),  (36),  (57)  and  (58),  we  have 


or 


-  -0  ■ 

or 

0 

8  " 

■  «■ 

6 

w 

0 

-  h  U 

0 

1 

0 

0 

*3 

r ' 9a 

L  *2 

-sp 

0 

sa 

1 

-8 

0 

-  -  sar 


L  *2 


v  h 


a 

0 

-sp 


0 

1 

0 


8 

0 

sa 


6 

0 

1 


which  implies  that 


-0  ■  w  h  (a  •  6  +  0) 


(59) 


r* 


and 


~  -  sa  -  w  h  (-s0  •  6  +  sa) 


Substituting  (50)  and  (51)  into  (59)  and  (60)  gives 


(60) 


L* 


wh5  Xl 

- - 7  TT 

1  ♦  wh(l  ♦  6Z)  a2 


and 


Ax'  ■ 


(61) 
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Ax'  ■ 


1  +  wh(l  ♦  62)  X2 


But  from  (51), 


i  5  A  *' 


Substituting  (61)  and  (62)  into  (63)  gives 


1  ♦  wh(l  ♦  52)  X2  1  +  vh(l  ♦  62  *2 


which  implies  that  X ^  -  -sX^.  Since  X^  and  X^  are  nonnegative  by  definition, 
we  have  -  X^«  But  this  contradicts  the  assumption  that  X^  i  X^.  There¬ 
fore,  the  necessity  part  is  proved.  Q.E.D. 

Theorem  III 

The  necessary  and  sufficient  condition  for  the  multiplicity  for  the 
singular  values  of  A  to  be  three,  i«e.,  X^  «  ■  Xj,  is  that  the  motion  consist 

of  rotation  around  an  axis  through  the  origin  only,  i.e.  Ax  ■  Ay  ■  Az  •  0. 

Also,  the  rotation  matrix  is  unique,  and  R  -  x/^A.  The  object  surface  can 
be  anywhere . 


Proof 


«  Xj,  then  (16)  gives 


x'2  +  y'2  +  *'2 
n  ■'n  n 


kM<xn  +  yn  +  *n> 


Since  any  3-D  rigid  body  motion  can  be  decomposed  into  rotation  followed 

by  translation,  we  first  rotate  the  object  such  that  (x  ,y  ,z  )  becomes 

n  n  n 

(xa,y4n,z") •  Then  we  carry  out  the  translation  which  changes  (x",  y",  z") 
—  —  —  n  n  n 


into  (x'»  y'»  z').  That  is, 
n  n  n 


+  Ay'  i 


where  R',  Ax',  Ay'  and  A*'  are  the  motion  parameters  in  the  (xq,  y^,  *n> 
space  as  defined  by  (23)  and  (24).  (65)  gives 


x«2  +  y**2  +  z"2 

a  'a  n 


y  zj  R*TF' 
l  n  'n  nj 


2  2  .  2 

x  ♦  y  +  z 
n  ^n  n 


This,  when  combined  with  (66),  gives 


*1  ♦  **  ♦  **  *  -  to’)2  ♦  CyJ  -  Ay')2  ♦  <»;  -  Az’)2 


From  (64)  and  (67),  we  have 


(k2\J-l)x^2  +  (k2\2-l)y^2  ♦  (k2\2-l)z^2  -  [2Ax’*x^  ♦  2Ay’*y]; 
+  2Az'*z^  -  Ax'2  -  Ay’2  -  As’2]  k2X2  -  0 


Since  (68)  is  true  for  all  x',  y '  and  s’,  by  equating  the  coefficients  of 
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all  powers  of  x',  y '  and  z'  to  zero,  we  have 
o  n  n 


Ay*  »  Az' 


1  or  k  *  t=- 

*1 


Therefore,  from  (24), 


Ax  ■  Ay  «  Az  «  0 . 


Then,  (9)  gives 


R  +  0  ■  k  A,  or  R  ■  X^  A 

Therefore,  we  ahve  proved  that  if  X^  «  X2  ■  Xj,  then  the  motion  consists 
of  rotation  around  an  axis  through  the  origin  only,  and  the  solution  for 
the  rotation  matrix  is  unique.  The  object  surface  can  be  anywhere.  This 
proves  the  necessity  part. 

We  now  proceed  to  prove  the  sufficiency  part.  If  the  motion  consists 
of  rotation  around  an  axis  through  the  origin  only,  i.e.,  Ax  •  Ay  »  Ax  ■  0, 


then  from  (9), 


A  -  k“l  R 


Let  U.  «  R,  V.  ■  I  and  A.  ■  k^I,  Then  (69)  becomes 


a  ■  uaaav; 


k 
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Sine*  and  are  orthonormal,  (70)  gives  the  SVD  of  A,  with  singular 
values  k  k  ^  and  k  Then  from  the  fact  that  the  singular  values  of 
any  matrix  are  unique,  we  see  that  A  has  three  identical  singular  values. 
This  proves  the  sufficiency  part.  Q.E.D. 

IV.  CONCLUSIONS 

Three  theorems  have  been  stated  and  proved  regarding  the  uniqueness 
and  the  computation  of  the  motion  parameters,  and  the  physical  descriptions 
and  classifications  of  the  actual  three-dimensional  motion  for  a  rigid  planar 
patch.  The  motion  parameters  are  unique  aside  from  a  scale  factor  for  the 
translation  parameters  if  the  singular  values  of  the  3x3  matrix  consisting 
of  the  8  pure  parameters  are  not  all  distinct,  otherwise,  the  number  of 
solutions  is  two.  The  distinction  between  the  cases  of  multiplicity  1  and  2 
lies  in  whether  or  not  the  translation  vector  coincides  with  the  normal 
direction  of  the  object  surface  at  If  there  is  no  translation  at  all, 

then  the  singular  values  are  all  Identical.  In  any  case,  once  the  eight 
pure  parameters  are  estimated,  which  can  be  done  by  solving  a  system  of 
linear  equations,  computing  the  singular  value  decomposition  of  a  3  x  3 
matrix  is  all  that  it  takes  to  obtain  the  3-D  motion  parameters  and  the 
directional  cosines  of  the  normal  direction  of  the  planar  patch. 
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Estimating  Three-Dimensional  Motion  Parameters 
of  a  Rigid  Planar  Patch,  III:  Finite  Point  Correspondences 
and  the  Three-View  Problem 
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ABSTRACT 


Two  results  are  presented  in  this  paper.  First,  it  is  shown 
that  in  estimating  three-dimensional  motion  of  a  rigid  planar  patch,  the 
eight  pure  parameters  used  in  [1]  and  [2]  are  uniquely  determined  from 
the  Image  correspondences  of  four  points,  no  three  collnear,  and  can  be 
estimated  by  solving  a  set  of  linear  equations.  The  second  result  con¬ 
cerns  the  three-view  problem.  It  is  proved  that  given  four  Image  point 
correspondences  in  three  perspective  views  of  a  planar  patch  undergoing 
general  three-dimensional  rigid  body  motion,  the  number  of  solutions  for 
the  motion  parameters  is  one,  as  opposed  to  two  [2]  when  only  two  perspec¬ 
tive  views  are  given. 
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I.  Introduction 

The  interest  in  motion  estimation  using  image  sequences  has  been  growing 
rapidly  in  many  fields  of  research  in  the  past  few  years.  The  efforts  in  the 
70's  were  primarily  focused  upon  two-dimensional  motion  estimation  [11-18]. 
Recently,  attention  has  been  gradually  shifted  toward  three-dimensional  motion 
estimation  [1-10,  19-21].  The  difference  between  2-0  and  3-D  motion  estimation 
is  not  just  the  degree  of  difficulty  or  complexity  in  solving  the  motion  equations. 
The  issues  of  uniqueness,  the  minimum  information  required  to  ensure  uniqueness 
and  the  3-D  structure  interpretation,  which  are  not  present  in  the  study  of  2-D 
motion  estimation,  make  the  study  of  3-D  motion  estimation  more  challenging  and 
interesting.  Furthermore,  due  to  the  nonlinearity  and  the  increase  of  the 
number  of  unknowns  of  the  motion  equations  for  3-D  motion  estimation,  the 
development  of  more  clever  and  efficient  ways  of  solving  the  motion  equations 
becomes  also  extremely  important. 

For  the  case  of  estimating  3-D  motion  of  rigid  curved  surfaces,  [3]  presen¬ 
ted  an  efficient  algorithm  for  determining  the  motion  parameters  exactly  without 
having  to  solve  nonlinear  equations,  and  was  the  first  to  analyze  the  problem 
of  how  many  image  point  correspondences  are  required  to  ensure  the  uniqueness 
of  motion  parameters. 

For  the  case  of  3-D  motion  estimation  of  a  rigid  planar  patch,  a  brief 
introduction  ia  given  in  Sec.  II.  In  this  paper,  it  is  proved  that  the  eight 
pure  parameters  [1,2]  in  the  two-view  problem  are  unique  given  the  image 
correspondences  of  four  points  no  three  colinear,  and  can  be  estimated  by 
solving  a  set  of  linear  equations.  For  the  three-view  problem,  it  is  proved 
that  given  four  image  point  correspondences  in  three  (distinct)  perspective 
views,  the  solutions  for  the  motion  parameters  are  unique. 


II*  The  Eight  Pur  a  Parameter*  and  the  A  Matrix 

The  baaic  geometry  of  the  problem  considered  in  [1]  and  [2]  is 
repeated  here  in  Figure  1.  Consider  two  frames  at  time  t^  and  t2  (t^  <  c^) 
For  a  particular  point  P  on  the  object,  let 

(x,y,z)  -  object-space  coordinates  of  a  point  P  at  time  t^. 
(x',y',z')  •  object-space  coordinates  of  P  et  time  tj. 

(X,Y)  -  image-space  coordinates  of  P  at  t^. 

(X’,Y')  -  image-space  coordinates  of  P  at  t^. 

It  is  obvious  from  Figure  1  that 


where  the  focal  length  is  normalized  to  one  for  convenience,  it  was  shown 
in  [1]  and  [2]  that  for  a  rigid  planar  patch  undergoing  3-D  motion  (a  rota¬ 
tion  with  an  angle  0  around  an  axis  through  the  origin  with  directional 
cosines  n^,  a,Jf  n^,  followed  by  a  translation  with  translation  vector  (Ax, 
Ay^z)) 


where  R  is  a  3  x  3  orthonormal  matrix  of  the  first  kind 


Fig.  1  Basic  geometry  for  three- 
dimensional  owe  ion 
eseiaecion. 


re-70»« 


n^  +(l-n^  )cos0 


(l-cosOJ+n^sinO 


R 


njn2 (l-co»9)-n3«in8 


2  2 

(l-co*0)+n3*ln0  n3  +  (1-02  )coa0 


(l-cos9)-nj*in0 


2  2 

n^n^l-cosOJ-^sin©  (l-cosOJ-m^ainO  n3  +(l-n3  )cos0 


Che  image-space  coordinates  before  and  after  the  notion  are  related  by 
a,X  +  a.Y  +  a, 


X' 


a^X  +  a„Y  +  1 


8 


Y'  - 


V  +  v  +  *6 


*7*  +  .SY  +  1 


where  the  a^'s  are  such  that  if  we  define  R  to  be  such  that 


2  2 

k  *  n3  +  (l-n3  )cos0  +  c.AZ 


and  let  ax  +  by  +  cz  »  1  be  the  equation  describing  the  object  surface 
before  motion  (at  t^),  then  ^ 

[a  b  c] 

► 

V  - 

The  eight  a^'s  in  A  are  called  the  pure  parameters  [1].  Given  A,  the 
actual  motion  parameters  can  be  obtained  by  simply  computing  the  SVD  of 
the  matrix  A  and  the  number  of  solutions  for  the  motion  parameters  is 
either  one  or  two  depending  on  the  multiplicity  of  singular  values  of  A  [2] 
As  for  the  uniqueness  of  the  pure  parameters  given  the  image  motion,  it 
was  shown  either  using  Lie  Group  Theory  [1]  or  elementary  algebra  [20] 
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Chat  given  tha  image  point  correspondences  of  the  vhole  plane,  the  pure 
parameters  are  unique.  It  is  the  purpose  of  Section  III  to  show  that 
it  takas  only  four  image  point  correspondences,  no  three  colinear  to  deter¬ 
mine  the  pure  parameters  uniquely.  The  algebraic  proof  in  Sectlonlll.  1 
is  direct  in  the  sense  that  the  coefficient  matrix  of  the  set  of  linear 
equations  that  the  eight  pure  parameters  must  satisfy  given  four  image 
point  correspondences  is  proved  to  be  nonsingular,  which  directly  leads 
to  the  conclusion  on  the  uniqueness  of  the  pure  parameters. -  The  geometrical 
proof  in  Section  III. 2  is  indirect  in  the  sense  that  the  pure  parameters 
are  not  directly  shown  to  be  unique,  but  rather  the  image  point  correspon¬ 
dences  of  the  whole  Image  plane  are  proved  to  be  fixed  given  four  Image 
point  correspondences.  In  order  to  ultimately  prove  that  the  eight  pure 
parameters  are  unique,  two  more  results  are  needed.  [1]  and  [20]  proved 
that  given  the  Image  point  correspondences  of  the  whole  plane,  the  eight 
pure  parameters  are  unique  if  the  3  x  3  A  matrix  is  nonsingular.  This  is 
the  first  result  needed.  The  second  one  is  Lemma  II  for  the  algebraic  proof 
in  Section  III. 1.  With  these  two  facts  and  the  geometrical  proof  in  Section 
III. 2,  one  can  conclude  that  the  eight  pure  parameters  are  unique  given 
four  image  point  correspondences,  noncolinear  both  before  and  after  the 
motion. 


In  order  to  prove  that  given  four  image  point  correspondences ,  the  pure 


parameters  are  unique,  we  take  the  following  approach.  Let  (X^'tY^')  be 
transformed  from  (X^,Y^),  i  ■  0,  ...»  3,  with  some  reference  pure  para* 


■r- 


X  ' 
*i 


*1  _  *2  Yj  3 


a  X  +  a  ^°^Y  +1 
*7  xi  +  *8  Yi  +  A 


V 


a  (o)x  +  a  (o)v  +  a  (o) 
*4  Xi  +  *5  Yi  +  *6 


»7(°>Xi  +  *8(°)y1  +  1 


Then  the  elements  of  the  matrix  M  in  (6)  contain  the  image  coordinates 


before  motion  (l.e.,  X^s '  and  Y^s 1 )  and  the  reference  pure  paramtera  a!  ’s 
Let  Ag  be  defined  as 


Ao  * 


(O) 

•a(0> 

V0) 

(o) 

•5(0> 

•6<8) 

(o) 

1 

It  is  to  be  shown  that  if  Ao  containing  the  reference  pure  parameters  is 

nonsingular,  then  give  four  image  point  correspondences  with  no  three  co- 

1 inear,  (5)  yields  only  one  set  of  solutions  for  the  a^ 's,  namely  the  refer 

<o). 

ence  pure  parameters  a{  s. 

In  order  to  simplify  the  analysis,  let  the  origin  of  the  Image 
coordinate  system  at  t^  be  located  at  one  of  the  image  points,  say  (Xq,Yq), 
at  t..  It  will  be  shown  at  the  end  of  this  section  why  this  simplification 


Substituting  (9),  (10)  and  (11)  into  (5)  givss 


(°).  (O)  .  -  (°) 

•o  •#! 


3  “8 


a  <°>-  (°)  -  «  (°> 
6  *7  6 

-  (o)  (o)  (o) 

6  *8  *5 


Note  that  B  doea  not  contain  the  unknown*  a^'*,  i  -  1,2, 4, 5, 7, 8.  Therefore 
if  D  la  nonalngular,  the  aolutlon  for  the  unknown* ,  *^»*2»*4»a5»a7»*a* 
unique.  Lean*  I,  below,  give*  the  exact  condition*  for  D  to  be  singular. 


Let  D  be  given  by  (13).  Then 

det(D)  -  djd^  det  (Ao)  det  (H)  (X^-X^)  (X^-X^) 


where 

Ao  Is  defined  in  (8) 

H  is  defined  in  (11) 

dl*d2*d3  are  given  in  (10) 

[Proof]  The  cofector  of  the  (ll)th  element  of  (13)  is  given  by 


cof(d1X1)  - 


After  some  straight forwerd  derivations,  the  above  becomes 

cof (djXj)  -  d2d3<X2Y3‘X3Y2)I,2(d3fl*dlf3)Y3'*3(d2fl“dlf2)Y21 
Similarly, 


d2Y2 

0 

0 

*2X2 

e2Y: 

d3Y3 

0 

0 

•A 

V: 

0 

dlxl 

dlYl 

fixi 

flY 

0 

d2X2 

d2Y2 

*2X2 

V: 

0 

djXj 

d3Y3 

f3X3 

V: 

and 


coftfjXj)  -  d143  0ElYj-X3Y1)[«1(d3f2-d2f3)Y3-«3(dlfJ-d2fl)Vl] 
cof(d3X3)  -  d2diat2Yl-XlY2)[.2(dlf3-d3fi)Yl-.1(d2(3-d3i2)Y2] 


Therefore, 


det(D)  -  d1X1cof(d1X1)-d2X2cof(d2X2)+d3X3cof(d3X3) 

•  dld2d3[,3(dlVd2fl>  <V3-W  (X3Yl'XlY3> 

•  *2  (dlf3*d3f l>  O'lVW  (X2Y3'X3Y2)<X3Yl'XlY3) 

+  VSVSV  (X1Y2-X2Y1)  (XjYj-XjYj)  (XjYj-XjYj) 

•  did2d3tel(d2f3“d3£2)“e2(dlf3*d3fl)+e3(dlf2"d2fl)1 

<X1Y2"X2Y1)  <X2Y3_X3Y2)  ^3Yl'XlY3> 

-  d1d2d3.det  (Aq)  (X1Y2-X2Y1)  (X2Y3-X3Y2)  (X^-X^) 

O.E.D. 

It  is  obvious  from  Lemma  I  that  if  d^,  d2  and  d3  are  never  zero 
(to  be  shown  later),  the  D  is  singular  and  only  if 

(i)  det(Ao)  •  0  (i.e.,  Aq  is  singular) 

(ii)  det (H)  *  0  (i.e.,  point  1,  2  and  3  are  colinear) 

(iii)  *iY2*X2Yl  *  ®  (!•*•»  point  0,  1  and  2  are  colinear) 

(iv)  ^2^3“^3^2  *  ®  (!•••»  point  0,  2  and  3  are  colinear) 

(v)  X3Yi*xiY3  “  0  (!•••»  point  0,  l  and  3  are  colinear) 


Note  that  (ii),  (iii),  (iv)  and  (v)  exhaust  ail  the  possibilities  for  any 
three  among  the  four  points  to  be  colinear.  We  now  show  that  d^,  d2  and  d 
are  strictly  positive. 


It  was  shown  in  [2]  that 


Substituting  (o)  into  the  above  gives 


zi'  “  k0^a7  ^ziXi  +  *8^2iYi  + 


a  (o)X  +  a  +  1 

*7  Xi  *8  *i  +  1 


—  k  ,  i 
2i  ° 


1,  2,  3 


(15) 


Since  the  object  points  must  be  in  front  of  the  camera,  z^'  >  1  (the  normalized 
focal  length)  and  z^  >1.  [2]  shows  that  kQ  <  0  corresponds  to  the  case  that 

the  object  points  move  to  the  back  of  the  camera.  kQ  obviously  cannot  be  zero, 
otherwise  (15)  would  imply  that  all  the  object  points  move  to  the  origin, 
therefore,  from  (16), 

2  ' 

a  +  a  +  1  ■  It  >0  ial,2  3. 

7  Ai  8  i  z^  o 


From  (10), 


Thus 
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di  >  0  for  i  ■  1,  2,  3. 

We  have  proved  Chat  D  in  (12)  Is  nonsingular  if,  and  only  if,  Aq  is  nonsingular 

and  no  three  points  among  the  four  are  co linear.  It  is  shown  in  the  following 

Lemma  that  the  restriction  on  A  being  nonsingular  need  not  be  imposed  If 

o 

none  of  the  three  points  are  colinear  both  before  and  after  the  motion. 

Lemma  II: 

Given  the  fact  that  the  image  points  before  motion  (at  t^)  are  not 
colinear  (or  equivalently,  at  least  three  points  are  noncolinear) ,  then  the 
following  three  statements  are  equivalent: 

(i)  kQ  is  sungular 

(11)  The  object  surface  passes  through  the  origin  at  t^. 

(ill)  All  the  image  points  after  motion  (at  t»)  are  colinear. 

[Proof]  We  prove  that  (i)  iff  (ii)  and  (ii)  iff  (iii). 

[(i)  (ii)] 


Let  the  SVD  of  A  be  given  by 
o 


where  X,,  X.  and  X.  are  the  singular  values  of  A  .  If  A  is  singular,  then 
1  Z  0  0 

one  of  the  singular  values  must  be  zero,  since  from  (17)  we  have 
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Vj  *  the  1th  column  of  V,j  »  1,2,3.  Since  U*  is  orthonormal,  (19) 
can  be  regarded  as  a  rotation  of  the  object-space  coordinate  system  around 
an  axis  through  the  origin  at  time  t^.  irom  the  fact  established  earlier 
that  one  of  the  X.s’  must  be  zero  if  A  is  singular,  we  see  from  (18)  that 

1  O  ' 

by  rotating  the  coordinate  system  around  an  axis  through  the  origin,  the 

object  surface  coincides  either  with  the  x"y"  plane  or  x"z"  plane  or  y"z" 

plane.  This  implies  that  before  rotating  the  coordinate  system  using  (19), 

the  object  surface  must  be  passing  through  the  origin. 

[(11) - )  (i)j 


Since  the  object  surface  passes  through  the  origin  at  t2»  by 
assumption,  (18)  becomes 


TT  Tt 

X1  V1J 
X2  v2  J 

T 

LX3  V3J 


(20) 


at  the  origin.  Were  Aq  to  be  nonsingular,  X^,  X^  and  X^  would  be  nonzero. 
Then  (20)  would  give 


mmm 

tmm 

0 

V 

T 

1 

0 

m 

V 

T 

2 

J  -  vTJ  or 

0 

V 

T 

3 

Lhs 

m  — 

ri 

p* 

0 

0 

X 

0 

0 

m 

0 

- 

J  or  J  ■ 

y 

- 

0 

0 

0 

z 

0 

m  mi 

L  — 

which  implies  Chet  Che  object  surface  at  t^  passes  through  the  origin. 

The  there  exists  a,  g  e  R  such  that 

z  »  a  x  +  g  y 

From  Equation  (0)  and  the  above, 
z-axz+gYz 

Thus 

a  x  +  g  Y  -  1 

which  implies  that  all  the  image  points  are  collnear  at  t^,  contradicting 
the  premise  of  the  lesma.  Thus  Aq  has  to  be  singular. 

[(iii)=*(ii)] 

From  (lii),  there  exists  er,  g  — *■  R  such  that 
Y'  -  a  X'  +  g 

From  Equation  (0)  and  the  above. 


or  or  x'  -  y'  +  g  z'  ■  0f  which  implies  (ii) 
l(ii)— »(iii)] 

From  (ii),  there  exists  or,  g  e  R  such  that 
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z-arx  +  gy-arxz+gYz 
or 

a  X  +  e  Y  -  1 
which  Implies  (111) . 


*  End  of  Proof  for  Lenina  II  * 

Up  Co  this  point,  we  have  proved  that  given  four  image  point 
correspondences  no  three  points  colinear  both  before  and  after  the  motion, 
the  pure  parameters  are  unique,  if  the  origin  of  the  image  coordinate 
system  coincides  with  one  of  the  image  points  at  t^.  It  is  shown  now 
that  the  latter  assumption  does  not  cause  any  loss  of  generality. 

First,  we  shall  show  that  given  four  image  point  correspondences 
in  two  frames,  one  can  derive  from  this  the  image  correspondences  of  the 
same  four  object  points  In  two  frames  taken  by  the  camera  in  the  same  posi¬ 
tion,  but  oriented  differently  such  that  one  of  the  object  points  coincides 
with  the  optical  axis,  i.e.,  the  z  axis.  Since  it  has  been  shown  earlier 
that  the  pure  parameters  are  unique  given  four  image  point  correspondences 
with  one  of  the  image  point  at  frame  1  at  the  origin,  we  see  that  the  3x3 
matrix  containing  the  pure  parameters  for  this  new  configuration  designated 
as  Ag,  is  unique.  Next,  we  shall  show  that  the  A  matrix  for  the  original 
configuration  is  similar  to  An,  and  can  be  determined  uniquely  from  Aq. 

The  proof  would  then  be  completed.  Now  we  furnish  the  details. 

Since  rotating  the  camera  is  equivalent  to  rotating  the  object 
points,  we  now  look  for  a  rotation  matrix  Rq,  which  can  rotate  the  point 
(xQ,  y0,  *0)  to  the  z  axis,  i.e.. 


Sine*  R  ,  if  norma lizfd,  wa  have 


I 


where  •  ■  (Xq  +  Yq  +  1)  and  Rq1,  R^  «f«  two  arbitrary  column  vectors 
such  that  Rq  is  orthonormal.  Note  that  although  Rq  is  not  unique,  any 
arbitrary  choice  of  Rq  will  lead  to  the  desired  conclusion,  as  to  be  seen 


later.  - 

It  is  to  be  shown  that  the  image  coordinates  of  the  four  points 

at  t.  and  t_  for  the  new  configuration  (all  the  points  are  rotated  by  R  ) 

can  be  derived  from  the  image  coordinates  for  the  original  configuration. 

Let  (xa£t  yni,  *ni)  be  the  object  coordinates  of  the  ith  point 

after  being  rotated  with  R  and  (x  .,  y  .)  be  its  image  corrdinates.  Then 

o  ni  ni 


R  T  J 
ol 

T 

*02  J 


•IW1  J 


where . J 


X  - 

ni 


T 

Rol  J 

•  t  Vol1 


T  -l 
Rol  J  «i 

*1* y.i]  J 


r  T  j' 
Rol  J 


i  [XqY^I]  j' 


Similarly, 


and 


R  T  j. 

Ynl- 


a[X  Y  1]  J* 
o  o 


X 

ni 


R  ^  J" 
Rol  J 


s[X  *Y  '  1]  J" 
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Y  «' 
ni 


Ro2T  J" 


«rxo'to,i]  j" 


where 


J"  - 


V 

V 


(24] 


(25] 


(26) 


From  (23)  -  (26),  we  see  that  2^,  Y^,  Xnll,  Y^’  are  functions  of  X±,  Yt, 

Xi'*  Xi*  Therefore,  the  image  point  correspondences  for  the  new  con¬ 

figuration  can  be  determined  directly  from  those  for  the  original  configuration. 


Since  1C  wee  proved  earlier  In  Chic  section  that  the  pure  parameters  are 
unique  given  four  Image  point  correspondences  with  one  of  the  point  at  the 


origin,  Aq  Is  unique.  Therefore,  comparing  (27)  and  (29), 


A  R 
o 


or 


A  -  R  A 
o  n 


(3 


Although  Rq  In  (22)  la  not  unique,  A  Is  still  unique  since  for  any  arbitrary 

choice  of  R  In  (22)  (Mote  that  different  R  would  result  in  a  different  A  ) 
o  on 

(30)  is  the  necessary  condition  for  all  possible  a'®.  Therefore,  the  pure 
parameters  for  the  original  configuration  are  unique. 

We  have  proved  that  given  four  image  point  correspondences,  no  three 
image  points  collnear  both  before  and  after  the  motion,  the  pure  parameters 
are  unique. 
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III. 2  Geometrical  Proof  for  the  Uniqueness  of  the  Point  Correspondences 

of  the  Whole  Image  Plane  Given  Four  Point  Correspondences 

It  is  worth  noting  that  the  geometrical  proof  presented  in  this 
section  does  not  lead  to  the  conclusion  on  the  uniqueness  of  the  pure  para¬ 
meters  directly.  As  explained  at  the  end  of  Section  II,  it  takes  Leona  II 
in  Section  III.l  and  the  results  in  [11(20]  Co  complete  the  proof. 

It  is  to  be  proved  that  given  the  correspondences  of  four  image 
points  In  two  perspective  views  no  three  colinear,  the  image  correspondence 
of  any  other  point  can  be  determined  uniquely.  In  particular,  let  A,  B, 

C,  D  and  E  be  five  arbitrary  points  in  frame  1,  such  that  no  three  are 
colinear,  and  let  A',  B',  C',  D*  be  the  given  corresponding  points  of  A, 

B,  C,  D  in  frame  2,  as  depicted  in  Figure  2.  We  would  like  to  show  that 
the  corresponding  point  E*  of  E  in  frame  2  is  uniquely  determined. 


A  .B 

•A' 

• 

•  B' 

.E 

C'# 

C  •£' 

•  • 

• 

C  D 

O' 

Frame  l(at  tt)  Frame  2  (at  t2) 


Figure  2  Five  point  correspondences  in  two  perspective  views. 

It  is  easy  to  see  that  if  a  set  of  points  in  the  image  space  are 
colinear,  the  corresponding  points  on  the  planar  patch  in  the  object  space 


must  also  be  coliaear  and  vise  versa.  In  fact,  let  (X^iY^),  i  -  1,  .  .  .  , 
n,  be  the  image  coordinates  of  n  coliaear  points  on  the  image  plane  and 
^Xi*yi,zi^*  *  ■  1*  ...  ,  n,  be  the  corresponding  points  on  the  planar  patch. 
Then  there  exist  a,  b  e  R  such  that 

X±  -  aYt  +  b 

Substituting  (o)  into  the  above  gives 

xi  “  ayi  +  b2i  < 

which  indicates  that  the  object  points  are  on  a  plane  passing  through  the 
origin.  Since  Af  B,  C,  D  are  not  colinear,  from  Lemma  II  of  Section  III.l, 
the  object  surface  cannot  pass  through  the  origin.  Therefore,  the  object 
surface  must  lie  on  the  intersection  of  the  object  surface  and  the  plane 
described  by  (£1).  Thus,  the  object  points  must  be  colinear.  The  converse 
is  obviously  true  since  we  can  regard  the  object  surface  as  the  image  plane 
and  vise  versa  and  then  repeat  the  above  argument. 

Next,  it  is  to  be  shown  that  given  the  correspondences  of  three 
colinear  points,  the  correspondences  of  all  the  other  points  on  the  line 
passing  through  these  three  points  are  determined. 

Consider  an  arbitrary  4th  point  on  the  line  passing  through  the 
given  three  points  in  frame  1.  Since  it  was  shown  earlier  in  this  section 
that  the  points  in  the  object  space  corresponding  to  a  set  of  colineer 
points  in  the  image  plane  must  also  be  colinear,  we  can  see  that  the  two 
sets  of  four  points,  one  set  on  the  image  plane,  the  other  set  on  the  planar 
patch,  are  in  perspective  correspondence  by  definition  [22].  Therefore, 
the  cross  ratio  [22,23]  of  the  four  points  in  the  image  plane  is  the  same 


as  the  cross  ratio  of  the  four  points  on  the  planar  patch  at  t^,  which 
remains  unchanged  from  t^  to  t^  since  the  object  undergoes  rigid  body 
motion.  Similarly,  the  two  sets  of  points,  one  on  the  image  plane  and  the 
other  on  the  planar  patch  at  ,  also  have  the  same  cross  ratios.  There* 
fore,  the  cross  ratios  for  the  two  sets  of  four  points,  one  at  t^  and  the 
other  at  t^,  are  the  same.  Then  by  definition  [22],  these  two  sets  of 
image  points  are  in  projective  correspondence.  It  is  well  known  in 
projective  geometry  (p.  83  [23])  that  the  projective  correspondence  between 
two  lines  is  fully  established  when  three  pairs  of  corresponding  points 
are  given.  Therefore,  we  can  see  that  given  the  correspondences  of  three 
colinear  points,  the  correspondence  of  any  other  point  on  the  line  is 
determined. 

Since  there  always  exist  two  straight  lines  not  parallel  to  each 
other  such  that  one  line  pauses  through  two  points  among  the  given  four 
points  A,  B,  C,  D  and  the  second  line  passes  through  the  other  two  points, 
we  can  assume  without  losing  generality  that  the  line  passing  through 
points  A,D,  denoted  by  A^,  is  not  parallel  to  BC.  Obviously,  A'D^  is 
also  not  parallel  to  B 1 C ^  in  this  case.  Since  none  of  the  three  among 
A,B,C,D  are  colinear,  the  point  lying  on  the  intersection  of  AD  and  BC, 
denoted  by  G,  does  not  coincide  with  any  one  among  A,  B,  C  and  D.  Similarly, 
the  point  lying  on  the  intersection  of  A'D^  and  Vd\  denoted  by  G' ,  does 
not  coincide  with  any  one  among  A',  B',  C'  and  D'.  If  E  lies  on  either 
one  of  AD  and  BC,  say  AD,  the  corresponding  point  E1  of  E  is  fixed,  since 
the  correspondences  of  the  three  points  A,  G,  C,  which  are  colinear  with 
E,  are  fixed.  On  the  other  hand,  if  E  does  not  lie  on  either  AD  or  BC,  let 


and  Lj  be  two  lines  in  frame  1  not  parallel  to  any  of  AD  and  as 
depicted  in  Figure  3,  and  let  the  points  at  the  intersections  of 
and  B?  be  denoted  by  H  and  J,  respectively,  and  the  points  at  the  intersec 

f4 

tions  of  L^,  1*2  and  AD  be  denoted  by  I  and  K,  respectively. 


Frame  1 


Frame  Z 


»-7390 


Figure  3.  The  point  correspondence  of  E  can  be  determined  from  the  point 
correspondences  of  A,  B,  C  and  D. 


Since  the  correspondences  of  B,  C  and  C  are  fixed,  the  correspondences 
all  the  points  on  the  line  BC  can  be  uniquely  determined.  Therefore, 
the  correspondences  of  H  and  J,  denoted  by  H*  and  J'  respectively,  are 
fixed.  Similarly,  the  correspondences  of  I  and  K,  denoted  by  I'  and  K* 
respectively,  are  also  fixed.  Therefore,  the  corresponding  point  of  E  in 


frame  2,  denoted  by  E',  which  lies  on  the  intersection  of  i'Hr  end 

is  fixed.  We  have  thus  proved  that  given  the  correspondences  of  four  image 

points  with  none  of  the  three  colinear,  the  correspondence  of  any  other 

point  in  the  image  plane  can  be  uniquely  determined.  Therefore,  the  mapping 

(x,y)  — *  (x  ,y')  is  fixed  for  all  (x,y)  e  R  .  Since  the  Image  points  are 

not  colinear,  according  the  Lemma  II  in  Section  III.l,  the  matrix  A  is 

o 

nonsingular.  Then,  from  [1][20],  the  pure  parameters  are  unique. 


IV.  Uniqueness  of  the  Motion  Parameters  Given  Four  Point  Correspondences 
in  Three  Image  Frames 

Consider  three  distinct  image  frames,  taken  at  three  time  instances 
t^,  t^  and  t^  (t^  <  t^  <  t,) ,  of  a  rigid  planar  patch  undergoing  three- 
dimensional  motion.  It  was  proved  in  Section  III  that  given  four  image 
point  correspondences  in  two  image  frames,  the  pure  parameters  are  unique, 
and  from  [2],  given  the  pure  parameters,  the  number  of  solutions  for  the 
real  motion  parameters  is  two  in  general,  unless  the  A  matrix  in  (4)  has 
multiple  singular  values.  In  this  section,  it  is  proved  that  with  four 
point  correspondences  in  three  image  frames,  the  solution  for  the  motion 
parameters  is  unique. 

Let  be  the  3x3  matrix  containing  the  eight  pure  parameters 
for  the  motion  from  t^  to  t^,  where  i  ■  1,  2,  3  and  J  ■  1,  2,  3  and  let 
k^  be  the  associated  constant  k  as  used  in  (4).  Consider  a  particular 
point  F  on  an  object.  Let 

(x,y,z)  -  object-space  coordinates  of  P  at 
(x',y',z')  -  object-space  coordinates  Of  P  at  t^. 

(x",y",z")  ■  object-space  coordinates  of  P  at  t^. 

(X,Y)  ■  image  space  coordinates  of  P  at  t^* 

(X'# Y')  ■  image  space  coordinates  of  F  at  t^. 

(X",Y")  -  image  space  coordinates  of  P  at  t^. 

It  can  be  shown  that  A. ,  ■  A.T*  and  k . ,  *  k 


In  fact,  from  (15), 


Since  four  image  point  correspondence*  with  no  three  colinear  are  given, 
according  to  Learns  II  in  Section  III. 1,  A^j  is  nonsingular.  Therefore, 
(32)  gives 


Comparing  (33)  and  (34)  shows  that 


is  one  possibility.  Sii.ce  it  was  proved  in  Section  III.l  that  given  four 
image  point  correspondences,  the  matrix  A  is  unique,  this  must  be  the  only 
possibility.  We  are  now  to  verify  the  following  composition  rules  for  A,. 


where  1  •  l,  2,  3,  j  ■  1,  2,  3,  n  »  1,  2,  3  end  n  +  1,  n  f4  j 
From  (15) »  we  have 
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Substituting  (38)  into  (39)  give* 


~x"' 

‘  x'' 

y" 

-  ^23^]2*23*12 

y' 

z" 

■»  a 

z» 

But  by  definition, 


A13  "  A23A12 


Ip  ■  V  If 

K13  *23*12 

Is  on*  posslbllty.  Sines  is  unique  given  four  image  point  correspondences, 
this  oust  be  the  only  possiblity.  As  for  other  velues  of  i,  j,  n  in  (36)  end 
(37),  the  proof  is  the  seme  except  for  the  change  of  indices. 

Since  moving  the  object  is  equivalent  to  moving  the  camera  so  far 
as  the  image  point  correspondences  are  concerned,  the  situation  can  be 
depicted  in  Figure  4  where  0^,  0^  and  0^  are  the  three  focal  points  for 
the  three  image  frames  when  the  planar  patch  is  considered  to  be  stationery 
while  the  camera  is  moving,  for  the  purpose  of  showing  the  relationships 


Since  four  image  point  correspondences  ere  given,  end  A^  *r«  fixed. 
Therefore,  ell  A^'s,  for  i  -  1,  2,  3  end  J  •  l,  2,  3  ere  fixed.  Then, 
from  [2],  for  the  motion  from  t^  to  t^»  there  ere  two  sets  of  solutions 
for  the  motion  parameters  given  A£^  end  for  the  motion  from  t ^  to  t^,  there 
ere  also  two  sets  of  solutions  for  the  motion  parameters  given  Since 

these  two  motions,  one  from  t ^  to  t^  and  the  other  from  to  t^,  can  be 
completely  independent  in  general,  the  only  possibility  for  the  solution 
of  the  motion  parameters  to  be  unique  is  that  not  both  of  the  two  solutions 
for  the  orientations  (i.e.,  the  directional  cosines  of  the  normal  directions 
of  the  object  surface)  of  the  planar  patch  corresponding  to  the  two  solutions 
of  the  motion  parameters  for  the  motion  from  t 2  to  t^  coincide  with  those 
for  the  motion  from  t 2  to  t^.  This  is  to  be  proved  by  contradiction. 

Assume  that  there  are  indeed  two  solutions.  Let  the  SVD  of  be 


Uij  Aij  Vij 


(42) 


where 


The  approach  we  shall  take  is  outlined  below: 


(i)  Prove  that 


3-34 


'21 


'23 


±1 


±1 


♦1 


(II)  Prove  that  ■  g  A^  for  some  constant  g. 

(III)  Provo  that  A^  ■  £*  I.  Than  tha  thraa  singular  values  for 
tha  notion  from  t^  to  t^  ara  Identical,  which  lnplias  from 
[2]  that  tha  solution  of  tha  motion  parameters  for  the  motion 
from  t^  to  t^  Is  unique,  contradicting  to  the  assumption 
that  there  are  two  solutions. 


The  details  for  tha  above  three  steps  are  now  given: 


(1)  Since  the  order  of  can  be  rearranged 

by  parmutatlng  tha  columns  of  U^j  and  In  (42),  we  can 
always  assume  that 

xf*l>  »  x»*l>  *  x<2*l> 

and 

X  (2,3)  *  .(2,3)  x  (2,3) 

1  a  X3 

If  any  of  tha  equality  signs  in  (44)  holds,  l.e.,  has  multiple  singular 
values,  then  for  the  motion  from  tj  to  t^,  the  solution  for  the  motion 
parameters  and  tha  orientation  of  the  object  surface  at  t ^  are  unique 
according  to  Theorem  I  In  [2].  Than,  for  tha  motion  from  t^  to  tj,  the 


(44) 


(45) 


solution  for  the  notion  paraaatsrs  oust  slso  bs  unique,  since  were  this 
false,  the  must  have  distinct  singular  values  and  there  would  be  two 
solutions  for  the  orientation  of  the  object  surface  at  t2>  Similarly, 
if  any  of  the  equality  signs  in  (45)  holds,  the  solutions  of  the  notion 
parameters  are  unique  for  both  the  notions  fron  t^  to  and  to  t^. 

If  only  the  inequality  signs  in  (34)  and  (35)  hold,  then  from  Theorem  II 
in  [2],  the  two  solutions  for  the  directional  cosines  of  the  planar  patch 
for  the  notion  from  t2  to  t^  are 


— x  a  A 

dl  “  bx  -  w  V21  0 


— V  A  A 

d  ■  K  m  w  V  0 

2  °2  w  V21  U 


and  for  the  motion  fron  t2  to  t^. 
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and 


"*2*" 

-6' 

V  * 

V 

■  v*  V 

W  v23 

0 

x 

1 

where  v  and  v'  are  some  normalizing  constants 


(49) 


Since  the  two  solutions  for  the  directional  cosines  of  the  object  surface 
at  t2  for  the  motion  from  to  t^  are  assumed  to  be  the  same  as  those  for 
the  motion  from  t2  to  t^»  either 


Ut  be  defined  as  the  nth  column  of  V..  for  n  ■  1,  2,  3.  Then, 

n  *j 


from  (46)  and  (47), 


Since  the  norm  of  7  Is  unity,  we  have  from  the  above,  • 


d  +  d 
1  +  2 


-  ±  »f  -l> 


II  t  + 


,(2.D 


We  now  show  that  *  Is  given  by  the  normalized  outer  product  of  and  d 


•  « 
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-5 

dl  X  d2  "  w  Vn 
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x  wvn 
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1 

•  • 

1 

-  »2  [S  V<2,1>  +  v^2,1)l  X  [-6  V^2,l)  +  V<2*1>] 

-  „2  [-  s2  v'2-1)  X  v«-l>  +  «  vf-l>  X  vf>2> 

-6vf’l>xvf’l)  +  v3<2-1>x  vf-^l 


where  "  X  "  stands  for  vector  outer  product 


With  either  (50)  or  (51),  we  heve  from  (52)  end  (57), 

end  from  (55)  end  (58),  V2(2,1)  -  ±  V<2,3)  end  from  (56)  end  (59), 

vf*1)  ■  +  v^2,3^.  Thus,  we  heve  proved  thet 


+1 

V21  "  V23  ±l 

+1 

(ii)  Let  (xn,  yn,  *n),  (xn’,  y^,  »n')  end  (xn»,  yn",  *n")  be  the 
new  coordinete  systems  obtained  by  roteting  the  coordinete 
systems  (x,y,s)  et  tj,  (x'.y',*')  et  tj  end  (x", y",*")  et  t3 
respectively,  using  the  orthonormel  metrix  V91  es  follows: 


3-41 


s.  x 
1  n 


*2  yn 


V23  y 


s.  * 

3  a 


where  s,  ■  +1,  i  ■  1,  2,  3.  Similarly,  from  (60)  ead  (63),  we  have 


•,  x 


•2  V  “  V23  y' 


3  n 


(64)  aad  (65)  Indicate  that  the  aew  coordinate  systems  (s  ,  y  ,  a  )  in 

nan 

(61)  and  (x  ",  y  ",  z  ")  in  (63)  can  also  be  obtained,  except  for  the  signs 
n  n  11 

by  rotating  the  old  coordinate  systems  using  V23  instead  of  V2^.  Note 
that  in  (64)  and  (65),  (x,  y,  z)  is  the  coordinate  system  at  t^  and  (x",  y",  z") 
is  the  coordinate  system  at  t^»  while  Is  the  matrix  in  (42)  containing 
the  singular  vectors  for  the  matrix  A^  which  characterizes  the  motion  from 
t2  to  t^.  Similarly,  in  (61)  and  (62),  (x,  y,  z)  is  the  coordinate  system 
at  tg  and  (x*,  y* ,  z')  is  the  coordinate  system  at  t^,  while  is  the  matrix 

in  (42)  containing  the  singular  vectors  for  the  matrix  A^  which  characterizes 
the  motion  from  t2  to  t^.  According  to  [2],  if  the  original  coordinate 
systems  in  the  object  space  are  transformed  as  in  (61)  and  (62)  using  V2l 
for  the  motion  from  t2  to  t^,  there  is  a  rigid  circle  lying  on  the  inter¬ 
section  of  the  object  surface  and  the  ellipsoid 


<2,l),2 


<2.D,2  2  .  rX<2.1),2  2, 


di  -  Vi  c*r  v  $ 


(66) 
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at  tg»  while  at  t^„  this  rigid  circle  oust  lie  oa  Che  intereection  of 
the  object  surface  and  the  sphere 

dl  *  *n’2  +  yn’2  +  *n'2  <”> 

where  d^  is  sons  constant.  We  would  like  to  show  that  this  rigid  circle 
oust  be  concentric,  on  the  object  surface,  with  another  rigid  circle  that 
lies  on  the  intersection  of  the  object  surface  and  the  ellipaoid 


*23  CPf’3’)2  ‘n2  +  *2<2,3)’2  yn2  +  P^’3’)2  \‘l 


at  t2  and  on  the  intersection  of  the  object  surface  and  the  sphere 
2  ••  2  _2 

d2  +*n  (69) 

■» 

at  t3  for  eome 

Because  of  (60),  (61)  end  (64),  we  can  see  that  the  principle 
axes  of  the  two  ellipeoids  in  (66)  and  (68)  are  the  sane.  From  (46),  (47) 
and  (61),  the  solutions  for  the  directional  cosines  of  the  planar  patch  at 
t2  in  the  new  coordinate  system  for  the  motion  from  t^  to  are  given  by 


Similarly,  from  (48),  (49)  and  (64),  the  solutions  of  the  directional  cosines 
of  the  planar  patch  at  tj  in  the  new  coordinate  system  for  the  motion  from 
tj  to  t3  are  given  by 


From  Che  above  cwo  equations,  we  can  aae  Chat  in  the  new  coordinate  system, 

the  normal  directions  of  the  object  surface  at  t^  oust  be  perpendicular  to 

the  yQ  axis,  both  for  the  motion  from  t^  to  and  for  the  motion  from  t^ 

to  t3»  Since  the  principle  axes  of  the  two  ellipsoids  in  (66)  and  (68) 

coincide  with  the  x  ,  y  and  z  axes  and  since  the  normal  direction  of  the 

n  n  n 

planar  patch  is  perpendicular  to  the  yQ  axis  (or  equivalently,  the  planar 
patch  is  parallel  to  the  y^  axis)  we  see  that  the  centers  of  the  rigid 
circles  lying  on  the  intersections  of  the  planar  patch  and  the  ellipsoids 
either  in  (66)  or  (68)  sust  be  on  the  xq  plane.  Obviously,  for  a 
particular  planar  patch,  as  d^  increases,  the  dimension  of  the  ellipsoid 
in  (66)  also  Increases  and,  consequently,  the  center  of  the  rigid  circle 
that  lies  on  the  Intersection  of  the  planar  patch  and  the  ellipsoid  becomes 
closer  to  the  x^  axis.  In  the  limit,  as  d^  goes  to  Infinity,  the  center 
is  on  the  xq  axis.  Similarly,  as  d2  becomes  large,  the  eenter  of  the  rigid 
circle  on  the  intersection  of  the  plsnsr  patch  and  the  ellispold  in  (68) 
approaches  the  s&  axis.  On  the  other  hand,  as  d^  decreases,  the  rigid 
circle  lying  on  the  intersection  of  the  planar  patch  and  the  ellipsoid 
in  (66)  gradually  shrinks  to  a  point.  For  a  particular  planar  patch,  let 
the  distance  between  this  limiting  point  snd  the  xq  axis  be  P^.  Similarly, 
as  decreases,  the  rigid  circle  lying  on  the  intersection  of  the  planar 
patch  and  the  ellipsoid  in  (68)  also  shrinks  to  s  point.  Let  the  distance 
between  this  limiting  point  snd  the  s_  axis  be  P„.  Then  it  is  seen  that 
the  distances  between  the  s_  axis  snd  the  centers  of  the  collection  of 
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r . 

r* 

r..* 


P. 


circles  lying  on  th*  planer  patch  and  tha  allipaoida  in  (66)  for  son*  range 
of  vary  between  0  and  P^,  while  the  diatancea  between  the  aQ  axis  and 
the  centers  of  the  collection  of  circles  lying  on  the  planar  patch  and  the 
ellipsoids  in  (68)  for  some  range  of  d 2  vary  between  0  and  P^.  Let  P^  be 
such  that 

0  <  p3  <  ain(plt  p2). 

From  the  above,  it  is  obvious  that  there  exist  at  least  two  rigid  circles 
lying  on  the  intersections  of  the  planar  patch  and  the  ellipsoids  in  (66) 
and  (68),  respectively,  such  that  the  centers  coincide  with  each  other  and 
the  distances  between  the  center  end  the  z  axis  is  p,. 

Let  the  equation  describing  the  object  surface  at  t£  be  expressed  as 


aX+Py+yi«l 

for  some  a,  0,  y  €  r.  Note  that  at  this  point  we  do  not  know  whether 
or,  8  and  y  in  (70)  are  unique  or  not.  However,  any  choice  of  or,  8  and 
Y  will  lead  us  to  the  conclusion  we  are  expecting  for  step  (li)  as  to  be 
seen  in  the  following.  Subetltuting  (70)  into  (66)  gives 


(70) 


(2,1)1* 

3 


(71) 


Similarly,  substituting  (70)  into  (68)  gives 
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Sine*  »n  no  laager  mun  la  (71)  nod  (72),  eh*  two  corn.  In  eh.  x  y 
ylo.  deecrlhed  by  (71)  end  (72)  *r.  the  v.rtlc.1  projection.  of 

ClrCl“  *»  «“  <V  yn.  «„)  .P*«  into  th.  x„  yB  plea*  end. 
therefor*,  net  be  tvo  elUpee.  thee  ere  Identical  up  to  e  .e.U  fector. 

*y  equating  the  coefficient,  la  (71)  and  (72)  up  to  .  poeltlv*  proper- 
tionality  constant  h,  hi  ha vs 


a  lV2’l)12  h  f  (V2’3T1 

-2*M  •  £ r 2on~J  | 


W2AV 


Lat  g 


V-  -  £  v2,sv 

7  r^i2  h^(2’3)j2 
/*, 

■g—  hj  .  (Nota  that  d^,  d2  «ad  h  ara  all  positiva). 


Than  (74) 


bacooMS 


(2,1) 


±  g  X. 


(2,3) 


But  ainca  cha  singular  valuaa  ara  nonnagaeiva  by  dafinition,  «a  hava 


Similarly,  (75)  gives 


*,«•«-  g  x2<2,3) 


Substituting  (78)  into  (73)  givss 


or 


or 


?*j<2,l)J2  ♦  [J*  Xj(2,3)]2  •  [g*j<2,3> 
C^2’15]2  -  [gX^2'3’]2 

H<2>1)  ■-  gx,«’s) 


From  (77),  (78)  and  (79),  vs  have 

iij,  -  g  a23 

(iii)  From  (36),  vs  have 


^3  "  A23  A12 


(35)  and  (81)  give 


13 


^3  1 


-  <CM  *23  T23)W2I  *2.  "gl*'1 
*  <023  *23  T2l><T21  *1}  "tf> 


Substituting  (60)  and  (80)  into  (82)  gives 
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U23  A23  V23  *  V23 


-1  -1  T 

«  *23  h  , 


«*1  D23  *33  ’  1 


A*1  n  T 
A23  U23 


-1  .  .-1  n  T 

g  a23  a23  u23 


(g*1  I)  0  * 


-  DAV' 


where 


A  A  g~l 


u  $  u23  ±1 


V  *  U„ 


Since  U23  is  orthonormal,  it  is  obvious  from  (85)  and  (86)  that  U  and  V 
are  orthonormal.  Also,  since  g  2  0,  A  in  (84)  is  a  diagonal  matrix  with 
nonnegative  diagonal  elements.  Therefore,  from  definition,  (83)  is  the 
singular  value  decomposition  of  A|3«  But  since  the  singular  values  are 


Since  the  three  eingular  values  of  for  the  notion  from  t^  to  ere 
identical,  fron  Theorem  III  in  [23,  the  motion  parameters  are  unique, 
contradicting  the  assumption  that  there  are  tvo  solutions.  We  have  thus 
proved  that  given  four  image  point  correspondences  in  three  image  frames, 
the  solution  for  the  motion  parameters  is  unique. 


V.  Conclusion* 


W«  have  shown  that  in  estimating  three-dimensional  motion  parameter* 
of  a  rigid  planar  patch  the  eight  pure  parameters  used  in  [1]  and  [2]  are 
unique,  and  can  he  determined  by  aolving  a  set  of  eight  nonsingular  linear 
equations  given  the  image  correspondences  of  four  points  with  no  three 
colinear  both  before  and  after  the  motion.  In  12]  it  was  shown  that  given 
the  eight  pure  parameters,  there  are  two  possible  solutions  to  the  motion 
parameters.  It  is  proved  in  this  paper  that  given  four  image  point  cor¬ 
respondences  in  three  (distinct)  perspective  views,  the  motion  parameters 
are  uniquely  determined. 
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